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Summary 

The  research  performed  under  the  contract,  during  the  period  17  March  1987  through  31 
March  1988,  can  be  divided  into  two  main  topics;  the  effects  of  ocean  continent  transition  zones  on 
L  -  g  waves,  and  modeling  regional  Love  waves  with  2-D  velocity  structures. 

In  section  II,  we  produced  hybrid  regional  and  telcseismic  SH  mode  sum  seismograms  after 
propagation  through  a  regional  transition  zone  or  other  heterogeneity  that  exists  as  part  of  a 
longer,  mostly  plane-layered,  path.  In  particular  the  effects  of  ocean  continent  transition  regions  of 
varying  lengths  on  the  propagation  of  Lg  waves  are  examined.  The  total  attenuation  of  Lg  due  tc 
geometry  for  propagation  along  a  partially  oceanic  path  is  at  most  a  factor  of  four.  Extending  the 
oceanic  path  length  changes  this  factor  to  at  most  6.  This  is  inadequate  to  explain  the  observed 
attenuation  of  Lg.  Thus,  additional  effects,  must  be  considered  to  provide  a  complete  explanation 
of  the  attenuation  of  Lg. 

In  section  III,  long  period  seismograms  recorded  at  Pasadena  of  earthquakes  occurring  along 
a  profile  to  Imperial  Valley  are  studied  in  terms  of  source  phenomena  versus  path  effects.  Some  of 
the  events  have  known  source  parameters,  determined  by  teleseismic  or  near-field  studies,  and  are 
used  as  master  events  in  a  forward  modeling  exercise  to  derive  the  Green’s  functions  that  describe 
the  propagation  effects  along  the  profile.  Both  timing  and  waveforms  of  records  are  matched  by 
synthetics  calculated  from  2-dimensional  velocity  models.  The  best  2-dimensional  section  begins 
at  Imperial  Valley  with  a  thin  crust  containing  the  basin  structure  and  thickens  towards  Pasadena. 
The  detailed  nature  of  the  transition  zone  at  the  base  of  the  crust  controls  the  early  arriving 
shorter  periods  while  the  edge  of  the  basin  controls  the  scattered  longer  period  surface  waves. 
Shallow  events  in  the  basin  are  easily  distinguished  from  deep  events  as  well  as  the  amount  of 
strike-slip  versus  dip-slip  motions  involved  from  the  waveform  characteristics  alone.  Those  events 
rupturing  the  sediments,  such  as  the  1979  Imperial  Valley  earthquake,  can  be  recognized  easily  by 
a  late  arriving  scattered  Love  wave  which  has  been  delayed  by  the  very  slow  path  across  the  shallow 
structure. 
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NUMERICAL  MODELLING  OF  SH  Lc  WAVES  IN  AND  NEAR 
CONTINENTAL  MARGINS 

Janice  Regan  and  David  G.  Harkrider 
Summary 

The  effect  of  transition  regions  between  continental  and  oceanic  structures  on  the  prop¬ 
agation  of  Ls  waves  from  continental  sources  is  examined.  In  particular,  the  attenuation  due 
to  variations  in  layer  thickness  in  such  transition  regions  is  calculated  and  explained  for  a 
suite  of  simple  models.  The  measured  attenuation  along  partially  oceanic  paths  beginning 
and  ending  in  a  continental  structure  is  at  most  a  factor  of  four.  This  is  inadequate  to 
explain  the  observed  extinction  of  L,  along  such  paths. 

The  method  used  to  calculate  the  results  presented  in  this  study  is  developed  and  its 
validity  and  accuracy  are  demonstrated.  Propagator  matrix  seismograms  are  coupled  into  a 
Finite  Element  calculation  to  produce  hybrid  teleseismic  SH  mode  sum  seismograms.  These 
hybrid  synthetics  can  be  determined  for  paths  including  any  regional  transition  zone  or  other 
heterogeneity  that  exists  as  part  of  a  longer,  mostly  plane-layered,  path. 

Numerical  results  presented  for  a  suite  of  continent  ocean  and  ocean  continent  transition 
models  show  several  trends.  Passage  through  a  continent  ocean  transition,  increases  the 
amplitude  and  coda  length  of  the  L,  wave  at  the  surface,  and  allows  much  of  the  modal 
Ls  energy  trapped  in  the  crust  to  escape  into  the  subcrustal  layers  as  body  waves.  The 
magnitude  of  both  these  effects  increases  as  the  length  of  the  transition  region  increases. 
When  the  wavefront  passes  through  the  oceanic  structure  following  the  continent  to  ocean 
transition  region  further  energy  escapes  from  the  crustal  layer,  and  produces  a  decrease  in 
Lf  amplitude  at  the  surface.  The  rate  of  amplitude  decrease  is  maximum  near  the  transition 
region  and  decreases  with  distance  from  it.  Passage  through  an  ocean  continent  transition 
region  causes  a  rapid  decrease  in  the  Lt  amplitude  at  the  surface  of  the  crust.  The  energy 
previously  trapped  in  the  oceanic  crustal  layer  spreads  throughout  the  thickening  crustal 
layer.  The  body  wave  phases  which  are  incident  on  the  continental  crust  in  the  ocean 
continent  transition  region  are  predominantly  transmitted  back  into  the  crust.  The  other 
body  wave  phases  reach  depths  below  the  depth  of  the  base  of  the  continental  crust  before 
reaching  the  ocean  continent  transition  and  thus,  escape  from  the  system. 

Introduction 

This  paper  presents  a  study  of  the  propagation  of  Lt  waves  across  ocean  continent  tran¬ 
sition  regions.  The  transition  regions  are  represented  by  simplified  models  each  consisting 
of  a  crustal  layer  with  a  smoothly  varying  thickness  above  a  half-space.  The  changes  to  a  I, 
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wavefield,  as  it  travels  across  such  a  transition  region,  are  modelled  using  a  hybrid  method 
which  combines  the  FE  method  and  the  Propagator  Matrix  technique  (PM)  (Harkrider. 
1964.  1970,  1981)  PM  seismograms  for  Lt  waves  from  a  continental  source  are  coupled  into 
a  FE  calculation  which  propagates  the  Lg  wavefield  across  the  continent  ocean  boundary. 
Results  from  the  Finite  Element  (FE)  calculation  are  then  coupled  into  a  second  FE  calcula¬ 
tion  which  propagates  the  Ls  wavefield  through  an  ocean  to  continent  transition  region.  The 
results  of  either  FE  calculation  may  be  propagated  through  a  region  of  horizontally  uniform 
waveguide  by  coupling  them  back  into  a  PM  calculation  using  the  Seismic  Representation 
Theorem  (RT)  (de  Hoop.  1958).  The  FE  to  PM  coupling  can  be  used  to  economical!}  in¬ 
vestigate  the  effects  of  long  ocean  path  lengths  between  regions  and  is  the  subject  of  later 
papers.  Here  we  restrict  ourselves  to  regions  in  and  near  transition  zones  separated  by  short 
(<  15(U-m)  ocean  paths. 

In  remainder  of  this  introduction  the  importance  of  this  work  to  major  areas  of  study  such 
as  tectonic  mapping  and  nuclear  discrimination  will  be  explained,  the  modal  interpretation 
of  Lg  on  which  the  calculation  of  synthetic  seismograms  is  based  will  be  justified,  and  the 
new  results  this  method  will  make  possible  will  be  discussed  and  related  to  previous  work. 

One  of  the  important  types  of  observational  studies  of  Lt  has  been  to  distinguish  regions 
with  oceanic  crustal  structures  from  those  with  continental  crustal  structures.  Press  and 
Ewing  (1952)  and  Bath  (1954)  observed  extinction  of  L,  when  the  propagation  path  included 
an  oceanic  portion  of  length  greater  than  200  km,  and  high  attenuation  or  extinction  when 
the  oceanic  path  length  was  as  short  as  100  km.  This  lead  to  the  commonly  used  approach, 
of  defining  paths  which  pass  under  oceans  but  do  not  attenuate  Lt  as  continental.  The 
results  of  the  present  paper  seriously  challenge  the  interpretation  that  paths  with  short 
oceanic  portions  which  show  little  or  no  L,  attenuation  are  necessarily  continental.  They 
may  necessitate  the  reassessment  of  some  of  the  results  of  studies  of  L,  in  many  regions  of  the 
world.  (Press  et  «J.,  1956;  Press,  1956;  Savarensky  and  Valdner,  1960;  Bolt,  1957;  Lehmann. 
1952,  1957;  Oliver  et  a/.,  1955;  Herrin  and  Minton,  1960;  Wetmiller,  1974;  Gregersen.  1984: 
Kennett  and  Mykkeltveit,  1984) 

Another  major  use  of  Lg  waves  is  the  determination  of  magnitudes,  of  explosions 

and  earthquakes.  Different  types  of  magnitudes,  including  are  compared  to  discriminate 
between  the  two  types  of  sources  (Blandford,  1982;  Pomeroy  et  ah,  1982).  7,  the  coefficient 
of  anelastic  attenuation,  which  is  derived  from  L,  wave  amplitudes,  is  important  in  studies 
of  m t,La  and  can  be  employed  to  assess  the  possible  destructiveness  of  earthquakes,  it  is 
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important  to  understand  if  reflections,  refractions,  or  diffractions  from  changes  in  crustal 
thickness  ,  generally  ignored  in  studies  measuring  m^L,  or  y  will  produce  significant  effect- 
not  accounted  for  in  the  interpretations  given  (Nuttli  1973,  1978,  19S1;  Herrmann  and  Xuttli. 
1982,  1975:  Street,  1976,  1984;  Street  et  al.,  1975;  Street  and  Turcotte,  1977;  Jones  et  ah, 
1977;  Bollinger,  1979;  Barker  et  ah,  1981;  Nicolas  et  ah,  1982;  Dwyer  et  ah,  19S3;  Chung  and 
Bernreuter,  1981:  Singh  and  Herrmann,  1983;  Campillo  et  ah,  19S4;  Herrmann  and  Kijko. 
1983). 

The  wavefields  transmitted  through  the  transition  region  models  are  calculated  to  model 
Ls.  The  preferred  interpretation  of  Lt  is  in  terms  of  a  superposition  of  higher  mode  surface 
waves.  The  higher  mode  surface  wave  interpretation  was  initially  unpopular  because,  in  its 
earliest  forms,  based  on  fundamental  mode  Love  waves  alone,  it  did  not  explain  the  vertical 
and  longitudinal  components  and  the  long  coda.  However,  it  subsequently  superceded  the 
alternate  explanation  in  terms  of  channel  waves  trapped  in  the  crust  above  a  low  velocity 
layer  for  the  following  reasons.  Oliver  and  Ewing  (1957,  1958),  Oliver  et  al  (1959).  and 
Kovach  and  Anderson  (1964)  showed  that  all  components  of  Lt  could  be  interpreted  b\ 
considering  both  higher  mode  Rayleigh  and  Love  waves.  Knopoff  et  al  (1975),  and  Panza 
and  Caicagmle  (i9/4,  i&75),  tfouchon  ( 1  S&u ,  iS$2)  used  the  higher  mode  interpretation  of 
L,  to  calculate  synthetic  seismograms  without  including  a  low  velocity  channel  below  the 
crustal  waveguide.  Other  phases  previously  defined  in  terms  of  the  channel  model  have  been 
successfully  modelled  using  the  higher  mode  surface  wave  model.  Schwab  et  al.  (1974),  and 
Mantovani  et  al  (1977),  considered  S0,  Panza  and  Caicagmle  (197o)  considered  Rg  and  L,. 
and  Stephens  and  Isacks  (1977)  considered  the  transverse  component  of  S„.  Clearly,  the 
multimode  surface  wave  explanation  of  L,  is  valid  and  useful.  However,  the  long  L,  coda 
observed  is  still  not  completely  understood  for  phase  velocities  less  than  2.8  km/sec.  The 
attribution  of  this  long  coda  to  diffraction  and  reflection  from  crustal  structure  is  supported 
by  the  results  presented  in  this  study. 

A  simple  parallel  of  the  multimode  surface  wave  interpretation,  which  is  a  very  useful 
aid  in  the  interpretation  of  the  wavefields  presented  in  this  study,  is  the  representation  of 
the  multimode  L}  arrivals  as  superpostions  of  multiply  reflected  post-critical  SH  and  SY 
rays  trapped  in  the  crustal  layer.  Bouchon  (1982)  used  this  type  of  interpretation  for  Ls 
arrivals  for  group  velocities  between  3.5  and  2.8  km/sec.  Pec  (1967)  and  Kennett  (19S6)  also 
used  the  r ay  approach  to  address  properties  of  Lt.  This  type  of  interpretation  can  also  be 
used  to  explain  where  structure  causes  conversions  from  Lt  to  distinct  body  waves  or  from 
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one  SJ1  mode  to  another.  Gregersen  (197S)  discusses  conversion  between  different  mo<ie_- 
of  Love  waves  and  between  Love  and  Rayleigh  waves  at  an  ocean  continent  boundaries 
Understanding  conversion  between  modes  of  L ,  and  between  Is  and  other  phases  is  an 
important  part  of  understanding  the  mechanisms  of  attenuation  of  I,  along  mixed  paths. 

Many  attempts  to  understand  the  propagation  of  seismic  disturbances  across  regions  of 
varying  structure  such  as  transition  zones  have  been  made.  First,  simple  models  were  used 
and  analytic  solutions  were  derived  for  soluble  special  cases,  then  increasingly  complicated 
models  were  considered  as  available  computational  power  increased.  The  types  of  models 
that  have  been  used  to  approximate  transition  regions  can  be  separated  into  several  types 
which  are  illustrated  in  Figure  1.  Sato  (1961a)  derives  analyse  results  for  models  of  type  la 
(Figure  la)  with  L=0  and  L>0.  Kennett  (1973)  has  considered  the  problem  of  seismic  waves 
interacting  with  a  layet  or  layers  in  which  propei ties  change  across  a  surface  perpendicular 
to  or  at  a  specified  angle  from  the  layering.  Several  types  of  solutions  for  models  of  type  lb 
(Figure  lb)  have  been  derived.  Sato  (1961b)  obtained  approximate  analytic  reflection  and 
transmission  coefficients,  for  h  <<  A,  the  incident  wavelength.  Hudson  and  Knopoff  (1964  ). 
KnopofF  and  Hudson  (1964),  Hudson  (1977),  and  Bose  (1975)  derived  similar  solutions  with¬ 
out  applying  the  h  <<  A  constraint.  Alsop  (1966)  developed  an  approximate  solution  for 
these  coefficients  applicable  when  all  energy  remains  in  Love  waves.  Gregersen  and  Alsop 
(1974,  1976)  extended  this  method  to  the  case  of  non-normally  incident  Love  waves.  They 
f^und  that  for  oblique  incidence  at  angles  less  than  forty  degrees  normal  incidence  is  a  good 
approximation.  Kazi  (1978a,  1978b)  derives  solutions  that  account  for  and  demonstrate 
the  importance  of  the  Love  waves  converted  to  scattered  body  waves  at  the  surface  step. 
Martel  (lybO)  used  a  FE  technique  and  spatial  filtering  to  isolate  the  diffracted  body  wave 
component.  Many  workers  have  studied  models  of  type  lc  (Figure  lc).  Knopoff  and  Mai 
(1967),  and  Knopoff  et  al  (1970)  derived  an  analytic  solution  usable  when  the  slope  of  the 
surface  (or  Moho)  in  the  transition  region  is  small.  Pec  (1967)  calculated  the  dispersion 
of  Love  waves  propagating  in  a  wedge  shaped  layer  and  found  that  the  largest  changes  in 
phase  velocity  and  amplitudes  at  short  periods.  Boore(1970)  studied  the  propagation  of  a 
simple  Love  wave  (T  >  20 sec),  modelled  as  a  Ricker  wavelet,  using  the  Finite  Difference  (FD) 
method.  He  noted  that  in  the  region  of  the  transition  that  mode  conversions  and  conversion 
to  other  types  of  waves  seemed  to  be  important.  Lysmer  and  Drake  (1971,  1972)  and  Drake 
(1972)  use  a  FE  method  based  on  Zienkiewicz  and  Cheung  (1967).  This  formulation  requires 
that  the  incident  modal  energy  is  exactly  equal  to  the  sum  of  the  reflected  and  transmitted 
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modal  energy  The  body  waves  present  in  the  system  produce  distortions  to  the  elastic  lavei 
over  a  half-space  eigenfunctions,  which  increase  for  higher  modes.  Drake  and  Bolt  (1980' 
used  the  same  method  with  a  model  of  type  Id  (Figure  Id)  to  study  a  more  realistic  model 
fitting  fundamental  mode  phase  velocity  data  for  events  normally  incident  on  the  California 
continental  margin  at  periods  between  4.4  and  60  seconds. 

All  the  studies  discussed  in  the  previous  paragraph  used  periods  much  longer  than  those 
that  will  be  considered  in  the  following  discussions.  The  shorter  periods  used  in  this  study 
allow  tin-  examination  of  the  effects  of  transition  regions  with  L  many  times  The  energy 
escaping  from  the  crustal  waveguide  is  shown  in  this  study  to  be  an  important  component 
of  the  explanation  of  the  attenuation  of  the  L,  phase  travelling  on  partially  oceanic  paths 
Previous  studies  considered  mainly  fundamental  mode  Love  wave  input  sampled  at  a  selection 
of  discrete  frequencies,  while  the  forcing  functions  used  in  this  study  are  a  sum  over  a  range 
of  frequencies  on  the  fundamental  and  first  five  higher  branches.  Of  the  previous  studies, 
only  Kennett  and  Mykkeltveit  (19S4)  have  generated  realistic  seismograms  similar  to  those 
used  in  this  study,  instead  other  studies  concentrated  on  measuring  phase  velocities  and 
transmission  and  reflection  coefficients. 

Propagator  to  Finite  Element  Coupling 

Within  a  layered  medium  the  "trapped”  wavefield  can  be  mathematically  constructed 
at  any  point  receiver  using  the  PM  technique  and  an  appropriate  form  of  a  source  represen¬ 
tation.  The  resulting  far-field  seismogram  will  include  not  only  a  direct  arrival  but  also  the 
superposition  of  many  multiple  critical  and  post-critical  reflections  which  produce  the  sur¬ 
face  waves  in  the  wavetrain.  To  simplify  the  graphical  representation  of  the  hybrid  method, 
shown  in  Figure  2,  each  such  seismogram  is  represented  by  a  single  solid  line  from  the  source 
to  the  receiver.  Coupling  a  wavefield  defined  by  the  discrete  sampling  of  such  seismograms, 
along  the  leftmost  edge  of  a  FE  grid,  into  that  FT  grid  is  straightforward.  The  seismograms, 
for  a  specified  source  function,  are  generated  at  a  group  of  receivers  equally  spaced  in  z 
located  a  horizontal  distance,  X,  from  the  source.  The  depth  spacing  between  the  receivers 
is  the  node  spacing  in  the  FE  grid  into  which  the  wavefield  is  to  be  coupled.  These  seismo¬ 
grams  are  applied  as  displacement  time  history  constraints  on  the  leftmost  column  of  nodes 
in  the  FE  grid,  thus,  completely  specifying  the  subsequent  motion  at  all  points  in  that  FF. 
grid. 

1  lie  seismograms  used  as  forcing  functions  can  be  generated  using  more  than  one  type 
of  sourff-.  For  a  line  source,  in  a  homogeneous  half-space  or  a  layered  half-space  the  applied 
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forcing  functions  are  uv(x,y,2).  For  a  double  couple  point  source  the  PM  solutions  can  be 
expressed  in  a  simple  form  when  only  variations  in  the  r  coordinate  are  being  discussed.  This 
form  is 


{v(r,<^,  2)}  =  v¥{4>,  2) 


dH?\kLr) 

dr 


where  v  =  2  for  a  strike  slip  source  and  v  -  1  for  dip-slip.  The  displacements  are  expressed 
in  cylindrical  coordinates  rather  than  in  the  Cartesian  coordinates  appropriate  to  the  R1 
integral  or  the  Cartesian  FE  representation.  When  these  displacement  seismograms  are  used 
as  input  forcing  functions  then  they  are  assumed  to  be  in  Cartesian  coordinates  (uv(x.  y. 
The  following  discussion  explains  why  the  hybrid  seismograms  produced  using  the  cylindrical 
forcing  functions  to  approximate  the  desired  cartesian  forcing  functions  approximate  the 
cylindrical  solution  {F(r,  <£,2)}  at  the  receiver. 

Consider  a  FE  grid  with  its  leftmost  edge  a  distance  X  =  r2  from  a  source,  and  a  receiver, 
where  hybrid  and  analytic  synthetic  results  are  recorded,  a  distance  XI  =  r2  from  the  source. 
Define  the  distance  propagated  within  the  FE  grid  as  Ar  =  r2  -  rj.  Since  we  are  considering 
the  case  of  r  large,  kL  »  *  {v(r2,^,x)}  can  be  expressed  in  terms  of  <*,  z,  n,  and  Ar  as 


Mr2,<£,2)}  =  {vfri.^.z)} 


cs  {v(ri,^,2)}  e_,ilAr  Ar  <<  ri  kirx  >>  1 
Now,  for  a  line  source  the  analogous  modal  continuation  relation  is, 


Comparing,  we  see  that  both  expressions  have  the  same  form.  In  each  case  the  displacement 
at  r2  can  be  expressed  as  the  displacement  at  n  multiplied  by  a  propagation  factor.  A  2  D 
FE  or  RT  calculation  will  give  the  same  propagation  factor  for  each  mode  as  the  analytic 
expression  above.  If  the  displacements  at  r2,  v(ri,4,z)  or  uv(z,z),  and  the  propagation  factor;. 
t->kL&r  or  are  correlated  by  considering  xj  =  rt  and  Ax  =  Ar,  then  F(r2,0,r).  the  3-D 

solution,  will  be  given  by  u„(x2,2),  the  result  of  the  2-D  FE  or  RT  calculation.  The  source 
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must  be  many  Ms,  kLr,  from  the  boundary  and  r  must  be  normal  to  the  FE  grid  edge  arid 
m  the  plane  of  the  FE  grid. 

The  PM  seismograms  fr(r,o,^)}  used  as  forcing  functions  in  the  tests  discussed  below 
show  the  overall  character  of  the  arrivals  seen  in  data  in  the  Lg  group  velocity  range  of 
between  2.8  and  3.5  km/sec  (Figure  3).  The  seismograms  are  complete  to  a  period  of  1.9 
sec  and  contain  most  of  the  energy  for  periods  down  to  0.5  sec.  The  heavy  rectangle  in 
Figure  3  labeled  input  window  delimits  the  time  window  used  to  select  the  portions  of  the 
seismogiams  used  as  displacement  time  history  forcing  functions.  The  forcing  function? 
used  in  subsequent  calculations  differ  from  the  seismogram  shown  in  Figure  3  in  that  no 
instrument  has  been  applied.  The  instrument  is  applied  to  the  FE  result. 

Designing  FE  Grids 

A  study  of  the  effects  of  the  length  of  simple  transition  regions  on  the  attenuation  of  SH 
type  L}  mode  sum  seismograms  passing  through  them  has  yielded  some  interesting  results. 
The  two  classes  of  transition  models  considered  are  illustrated  in  Figure  4.  Calculations  were 
performed  for  four  individual  models  from  each  class,  for  a  continental  reference  model,  and 
for  an  oceanic  reference  model.  The  difference  between  individual  transition  models  was. 
L,  the  length  of  the  transition  region.  Real  ocean  to  continent  type  transitions  occur  over 
lengths  of  50  to  300  km  (Keen  and  Hyndman,  1979;  Hinz  et  a/.,  1979;  Eittreim  and  Grantz. 
1979;  Le  Douaran  et  sJ.,  1984).  However,  an  upper  limit  on  L  of  100  km  was  imposed  by 
computation  limits.  Thus,  the  lengths  used  for  this  investigation  were  a  step  transition  (0 
km).  25,  50.  arid  100  km.  In  order  to  discuss  the  results  of  the  FE  calculations  using  these 
models  one  must  first  describe  the  models.  Also,  the  methods  used  to  obtain  and  display 
the  results  of  the  calculations  using  those  models  must  be  considered. 

The  first  class  of  models  are  used  to  describe  continent  to  ocean  transition  regions.  In 
further  discussions  these  models  will  be  referred  to  as  C/O  models,  and  the  transitions  they 
represent  as  C/O  transitions.  Similarly,  the  second  class  of  models,  used  to  describe  ocean 
to  continent  transition  regions,  will  be  referred  to  as  O/C  models  and  the  transitions  they 
represent  as  O/C  transitions.  The  continental  reference  model  consists  of  a  32  km  thick 
crustal  layer  over  a  half-space.  The  oceanic  reference  model  consists  of  two  5  km  layers,  one 
water  and  one  crust,  over  a  half-space.  Each  C/O  or  O/C  transition  region  is  characterized 
by  a  continuous  rate  of  thinning  or  thickening  of  the  crustal  layer.  In  all  models  the  crustal 
laver  has  a  SH  wave  velocity,  of  3.5  km/sec  and  a  density  of  2.7  g/cc,  while  the  half-space  has 
an  SH  wa’  <■  velocity  of  4.5  km/sec  and  a  density  of  3.4  g/cc.  Each  class  of  transition  model4 
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has  the  same  DC’s  applied  to  each  of  its  members.  The  forcing  functions  used  to  drive  the 
C/0  transition  and  continental  reference  FE  calculations  are  a  vertical  section  of  sixty  mode 
sunn  seismograms,  calculated  at  depth  intervals  of  0.5  km  beginning  at  the  surface,  for  a 
strike  slip  source  at  8  km  depth  at  a  distance  A'  =  1500  km.  The  forcing  functions  for  the 
O/C  transition  and  oceanic  reference  tests  are  recorded  during  the  50  km  C/O  transition 
calculation.  They  consist  of  a  depth  section  of  hybrid  seismograms  recorded  25  km  past  the 
oceanic  end  of  the  50  km  C/O  transition  region,  1755  km  from  the  source.  The  vertical 
spacing  within  the  depth  section  is  0.5  km.  All  other  BC’s  are  identical  for  all  models.  The 
transmitting  BC  is  applied  to  all  edge  nodes  of  the  grid,  excluding  only  the  free  surface  and 
the  nodes  to  which  forcing  functions  are  applied. 

The  first  step  in  designing  a  grid  for  FE  calculations  is  to  determine  grid  spacing,  dx. 
and  the  time  step  duration,  dt.  It  is  important  to  maximize  dx  and  dt.  At  least  six  nodes 
per  wavelength  are  needed  to  avoid  numerical  dispersion  problems.  To  maintain  numerical 
stability,  the  wavefront  can  travel  no  more  that  half  the  grid  spacing  per  time  step.  Thus, 
dx  and  dt  must  be  chosen  to  satisfy 


dx  < 


dt  < 


dx 


~  2VL 


where  f  is  the  highest  frequency  in  the  waveform  to  be  modelled  and  Vkiin  and  V(mar)  are 
the  minimum  and  maximum  S  velocities  respectively  in  the  dominant  part  of  the  model.  In 
this  study,  we  are  considering  L,  waves  with  a  predominant  period  of  approximately  1  sec. 
in  a  medium  with  Vmin  =  3.5  km/s  and  Vmax  =  4.5  km/s.  Thus,  we  have  chosen  dx=.5  km 
and  dt=.05  sec. 

The  next  step  in  designing  the  grid  is  determining  the  dimensions  of  the  grid,  nx  and 
ny,  the  duration  of  the  input  forcing  functions,  Tm,-,  the  duration  of  the  calculated  time 
series,  Tmcau,  and  the  location  of  the  transition  region  within  the  grid.  The  location  of  the 
transition  region  within  the  grid  is  defined  in  terms  of  the  distances  from  the  leftmost  grid 
edge  to  positions  A,  B,  C,  D,  E,  in  Figure  4.  The  values  of  these  parameters  were  chosen  to 
satisfy  two  criteria. 

(A)  a  seismogram  of  duration  D,  seconds  can  be  recorded  at  A  (Figure  4)  before  the  multiple 
reflection  of  the  input  wave  from  the  beginning  of  the  transition,  B,  to  the  leftmost  grid 
boundary  and  back  again  reaches  A. 
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iB)  a  seismoeram  uncontaminated  by  the  multiple  reflection  with  duration  D,  seconds  can 

be  recorded  at  the  receiver  closest  to  the  rightmost  edge  of  the  grid.  This  receiver  is 

defined  to  be  at  a  distance  xt  from  the  leftmost  grid  edge. 

Criteria  A  and  B  concern  themselves  only  with  reflections  from  the  leftmost  edge  of  the  grid. 
Nonphysical  reflections  from  the  bottom  and  the  rightmost  edge  of  the  grid  are  removed 
using  transparent  BC's  which  will  be  discussed  later.  The  duration  D,  was  chosen  to  be  55 
seconds,  the  observed  coda  length  for  a  SH  Ls  mode  sum  seismogram  at  a  distance  of  1000 
km  from  the  source.  For  all  the  models  used  in  this  study  the  values  of  parameters  defining 
the  size  of  the  grid  and  the  location  of  the  transition  region  within  it  are  given  in  table  1 
Distances  are  given  as  the  number  of  nodes  in  the  horizontal  direction  from  the  left  edge  of 
the  grid  to  the  depth  section  or  boundary  indicated.  The  model  names  indicate  L.  then  the 
model  type  (f=C/0,  r=0/C,  ref=reference). 

Sampling  FE  Solutions 

Analysis  of  the  effects  of  various  transitions  on  the  waveforms  and  amplitudes  of  L,  waves 
using  FE  techniques  requires  that  the  motions  of  the  nodes  of  the  FE  grid  be  sampled  so  that 
the  progress  of  the  L ,  waves  across  the  transition  can  be  observed.  Two  methods  of  sampling 
are  used  in  this  study.  Complete  displacement  time  histories  are  recorded  for  selected  nodes, 
and  the  displacements  of  all  nodes  in  the  grid  are  recorded  at  given  time  intervals.  The  first 
approach  produces  seismograms  which  can  be  used  to  illustrate  variations  of  amplitude  and 
waveform  with  distance  or  depth,  the  second  approach  produces  time  slices  and  is  a  clear 
way  to  illustrate  the  propagation  and  distortion  of  wavefronts  caused  by  passage  through 
the  inhomogeneous  structure. 

For  each  model  time  slices  were  recorded  every  250  or  400  time  steps.  Displacements 
are  graphically  represented  by  symbols  plotted  at  an  array  of  points  depicting  the  nodes 
in  the  FE  grid.  The  size  of  the  symbol  plotted  at  the  node  is  increased  as  the  absolute 
value  of  the  displacement  increases,  producing  darker  areas  where  larger  displacements  are 
occurring.  Each  of  the  time  slices  is  self  scaled,  that  is  the  largest  value  of  the  absolute 
value  of  amplitude  in  the  grid  is  set  to  the  largest  symbol  size.  At  all  other  nodes  the  ratio 
of  the  nodal  amplitude  to  the  maximum  amplitude  is  used  as  a  factor  to  scale  the  symbol 
size.  A  minimum  amplitude  cutoff  is  defined  below  which  no  symbol  is  plotted.  For  the  time 
slices  displayed  in  this  paper  the  cutoff  is  2%  of  the  maximum  amplitude.  Thus,  the  smallest 
symbols  cover  a  range  of  amplitudes  between  1%  and  10%  of  the  maximum  amplitude.  The 
self  scaling  of  the  time  slices  means  that  successive  time  slices  may  show  the  same  absolute 
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amplitude  as  a  different  symbol  size.  Thus,  the  same  region  of  the  waveform  will  appe?.; 
darker  on  a  time  slice  with  a  given  maximum  amplitude  than  on  another  time  slice  wit;, 
a  larger  maximum  amplitude.  This  difference  must  be  remembered  when  interpreting  the 
time  slices. 

For  each  model  one  surface  section  and  several  depth  sections  of  seismograms  were 
recorded.  Each  surface  section  contains  seimograms  recorded  at  intervals  of  5  km  along 
the  solid  surface.  Each  depth  section  contains  a  set  of  seismograms  recorded  in  the  same  col¬ 
umn  of  the  FE  grid,  starting  in  the  upper  row  of  the  grid  and  equally  spaced  in  depth  below 
it.  Depth  sections  with  a  vertical  spacing  of  2.5  km,  were  recorded  at  distances  including 
those  corresponding  to  positions  A  through  E  (Figure  4)  for  each  model.  For  the  each  refer¬ 
ence  model,  depth  sections  were  recorded  at  the  distances,  X2,  used  for  each  corresponding 
transition  model. 

The  amplitude  of  the  first  large  positive  and  negative  peaks  in  L,  mode  sum  seismograms 
(figure  3)  are  extremely  sensitive  to  distance  from  the  source,  and  to  the  time  spacing  and 
starting  time  chosen  when  calculating  the  synthetics.  It  is  not  unusual  for  very  small  changes 
in  these  parameters  to  produces  variations  in  amplitude  of  several  percent  .  An  improved 
filtering  technique  would  reduce  the  variation  slightly  by  removing  more  of  the  sampling 
effects  due  to  the  non-zero  spectral  amplitude  at  the  nyquist  frequency.  However,  a  significant 
portion  of  the  variation  is  apparently  due  to  changes  in  the  patterns  of  interference  between 
multiple  arrivals.  In  addition,  the  direct  S  pulse  and  the  end  of  the  S„  wavetrain  arrive  at 
the  same  time  as  the  initial  pulses  in  the  L,  wavetrain.  The  longer  period  portions  of  the 
S  and  S„  arrivals  have  been  removed  by  the  band  pass  filter.  However,  some  non  Ls  energy 
will  be  present  in  the  first  few  cycles  of  the  trace.  Unfortunately,  the  initial  high  frequency 
peaks,  which  are  effected  by  the  problems  discussed  above,  yield  the  maximum  peak  to 
peak  amplitude.  This  indicates  that  peak  to  peak  amplitude  will  be  a  poor  measure  of  the 
L,  amplitude.  A  more  stable  measure  of  amplitude  correspondence  between  seismograms 
of  this  type  is  the  RMS  amplitude  calculated  over  some  time  window  appropriate  to  the 
seismograms  being  compared. 


RMS  = 
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where  m  is  the  number  of  points  in  a  sampling  window  of  duration  Trms  seconds.  Such  an 
RMS  amplitude  measure  will  reduce  the  effects  of  S  or  Sn  contamination  and  of  any  other 
instabilities  effecting  the  initial  arrivals. 

The  evaluation  of  RMS  amplitudes  is  straightforward  once  the  sampling  window  has 
been  chosen.  The  location  of  the  sampling  window  with  respect  to  the  arrivals  of  maximum 
amplitude,  and  the  duration  of  the  trace  contained  within  the  sampling  window,  have  a 
significant  effect  on  the  value  of  the  RMS  amplitude.  Care  must  be  taken  to  choose  windows 
for  two  sets  of  results  that  produce  meaningful  comparisons.  When  results  of  two  separate 
FE  calculations  are  being  compared,  at  corresponding  nodes,  both  seismograms  begin  at 
the  same  absolute  time.  When  successive  seismograms  in  a  depth  or  distance  section,  from 
a  single  calculation,  are  being  compared,  the  duration  of  the  sampled  portion  of  the  trace 
with  negligible  amplitude  that  occurs  before  the  first  arrival  must  be  constant.  A  general 
approach  which  does  not  require  the  use  of  theoretical  travel  times  was  chosen  to  determine 
windows  in  the  latter  case.  Two  different  methods  of  selecting  the  first  and  last  point  in  the 
sampling  window  have  been  used.  The  first  method  selects  by  bracketing  the  portion  of  the 
seismogram  with  significant  amplitude  according  to  the  following  algorithm, 

(1)  Set  a  cutoff  value  for  the  amplitude  at  some  fraction  of  the  maximum  absolute  value  of 
amplitude  in  the  seismogram  (usually  .01  or  .05) 

(2)  Let  the  sample  at  the  location  of  the  maximum  absolute  value  be  the  first  sample  in  each 
of  two  series.  One  series  proceeds  forward  in  time,  the  other  backward.  Scan  each  series 
until  a  subseries  of  samples  n  seconds  in  duration,  all  with  amplitudes  smaller  than  the 
cutoff  amplitude,  is  found.  The  first  point  in  each  of  these  subseries  defines  an  endpoint 
of  the  longest  possible  sampling  window. 

(3)  Calculate  the  RMS  amplitude  for  sampling  windows  with  durations  15,  20,  25,  ..,  85  sec. 
If  the  duration  of  a  sampling  window  exceeds  the  duration  of  the  seismogram  then  the 
longest  sampling  window  used  is  the  duration  of  the  seismogram. 

This  algorithm  fails  when  a  period  of  sustained  low  amplitude  (>  nsec )  occurs  within  the 
desired  trace,  or  when  small  arrivals  precede  the  main  arrival  on  only  a  subset  of  the  traces. 
The  second  algorithm  selects  the  point  in  the  time  series  for  a  time  ten  seconds  before 
the  peak  maximum  absolute  amplitude  as  the  first  point  in  the  window.  The  end  of  the 
seismogram  is  defined  by  the  first  point  in  a  five  second  sample  of  sustained  low  amplitude. 
The  combination  of  these  two  methods  gives  stable  RMS  amplitudes  in  all  cases  considered. 


n 


These  algorithms  yield  a  series  of  RMS  amplitude  values  that  can  be  used  to  compare 
results  derived  using  different  methods.  If  RMS  amplitudes  agree  within  all  period  ranges, 
and  waveforms  are  similar  the  fit  is  considered  to  be  excellent.  Examining  the  behavior  of 
the  RMS  amplitude  as  a  function  of  sample  window  length  can  give  insight  into  the  nature 
of  and  the  underlying  reasons  for  misfit  between  methods. 

Tests  of  Accuracy:  SH  Lt  Mode  Sum  Input 

A  series  of  tests  were  conducted  to  verify  the  accuracy  of  the  coupling  method  and  the 
efficiency  of  the  absorbing  BC’s.  The  tests  used  the  continental  reference  model  and  a  set 
of  ninety  forcing  functions  calculated  for  a  strike  slip  point  double  couple  source  at  a  depth 
of  8  km  and  a  distance  X=1500  km.  Transparent  BC’s  were  applied  when  necessary,  as 
explained  in  the  discussion  of  the  test  results.  The  BC’s  are  discussed  below  before  the  test 
results  are  presented. 

The  BC’s  are  introduced  to  remove  the  nonphysical  reflections  created  by  the  artificial 
grid  boundaries  at  the  ends  and  bottom  of  the  grid.  Removing  these  reflections  by  extending 
the  grid  requires  an  increase  in  execution  time  of  between  30%  and  500%.  As  a  comparison 
the  approximate  BC  used  in  this  study  increases  execution  time  by  3-5%.  The  BC  used 
here,  suggested  by  Frazier,  Alexander,  and  Petersen  (1973),  calculates  the  average  of  the 
rigid  boundary  and  the  free  boundary  displacements  for  each  edge  node  at  each  time  step 
For  a  normally  incident  plane  wave  this  average  exactly  represents  the  transparent  boundary. 
However  in  practice  the  incident  wavefront  is  neither  normally  incident  nor  a  plane  wave. 
This  means  that  the  actual  value  at  the  transmitting  boundary  is  a  linear  combination  of 
the  rigid  boundary  and  free  boundary  solutions  whose  coefficients  depend  upon  the  angle  of 
incidence  of  the  energy.  The  boundary  condition  used  here  assumes  that  the  average  of  the 
two  solutions  will  in  most  cases  be  the  best  approximation  to  the  transparent  boundary  that 
can  be  simply  implemented. 

The  results  of  the  accuracy  and  BC  tests  presented  below  are  illustrated  in  Figure  5. 
A  continental  reference  structure  is  used  for  both  the  FE  and  PM  portions  of  the  path. 
This  allows  the  direct  comparison  of  the  hybrid  solutions,  those  solutions  propagated  from 
the  source  to  the  receiver  using  a  combination  of  methods,  to  direct  analytic  synthetics, 
those  calculated  entirely  with  the  PM  technique.  Four  separate  calculations  were  performed. 
Results,  from  each  of  these  calculations,  were  recorded  at  20  and  40  km  into  the  grid. 
First,  mode  sum  direct  analytic  seismograms  for  the  same  source  used  to  generate  the 
C/O  transition  forcing  functions  were  calculated.  Next,  seismograms  were  generated  using 
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the  hybrid  method  in  a  grid  long  enough  (100X50  km)  to  prevent  contamination  from  end 
reflections.  Third,  seismograms  were  generated  using  the  hybrid  method  and  a  short  (50X50 
km)  grid  with  the  transparent  BC  applied  at  each  node  on  it's  rightmost  edge.  Finally, 
the  previous  FE  calculation  was  repeated  without  the  transparent  BC’s.  The  excellent 
agreement  between  the  analytic  and  hybrid  synthetics  verifies  the  validity  of  the  coupling 
method  applied  to  a  layered  half-space.  Small  discrepancies  are  seen  in  the  higher  frequency 
component,  particularly  in  the  first  twenty  seconds  of  the  trace.  In  the  last  trace  in  each  group 
in  Figure  5  the  single  and  multiple  reflections  from  the  grid  edge  are  clear.  Comparison  of  the 
long  grid  and  the  short  grid  with  and  without  BC’s  shows  that  most  of  the  reflected  amplitude 
has  been  removed  by  the  BC's.  The  misfit  is  lower  for  the  longer  period  component  of  the 
traces.  The  significant  misfit  is  coincident  with  reflections,  and  is  largest  for  the  multiple 
reflections.  Despite  visible  differences  in  waveform  the  seismograms  shown  in  Figure  5  have 
RMS  amplitudes  that  agree  to  within  less  than  2%  for  all  RMS  window  lengths.  This 
indicates  that  small  changes  in  waveform  may  be  expected  but  the  RMS  amplitudes  of  the 
seismograms  should  be  stable  and  not  significantly  contaminated  by  incompletely  removed 
reflections  from  the  grid  edges.  The  increased  discrepancies  in  both  waveform  and  amplitude 
introduced  by  the  multiple  reflections  will  be  avoided  in  the  transition  FE  grids  described 
below.  This  reduces  the  discrepancies  in  RMS  amplitude  to  less  than  1%. 

Reflections  from  the  bottom  edge  of  the  grid  should  also  be  considered.  The  transparent 
BC  can  be  very  inefficient  for  the  small  angles  of  incidence  seen  at  the  grid  bottom  when 
a  distant  source  is  considered  (Regan,  1987).  Fortunately,  in  a  layered  half-space,  where  Lf 
waves  can  be  modelled  as  the  superposition  of  multiply  reflected  post-critical  SH  waves  in 
the  crust,  most  of  the  energy  in  the  SH  type  L,  waves  should  be  trapped  in  the  crustal  layer. 
The  time  slices  discussed  later  (Figures  6,  12)  show  that  the  SH  type  L,  waves  propagating 
in  a  layer  over  a  half-space  are,  in  the  most  part,  confined  within  the  layer.  RMS  amplitudes 
agree  to  within  less  than  0.1%  with  and  without  the  BC’s  (Regan  1987).  However,  the 
purpose  of  the  present  study  is  to  examine  the  effects  of  continental  oceanic  boundaries  on 
the  transmission  of  L,  mode  sum  seismograms.  When  the  crustal  layer  is  thinned  or  thickened 
with  distance,  the  modes  are  no  longer  completely  trapped  within  the  layer.  Energy  can  be 
converted  to  modes  compatible  with  the  local  layer  thickness  and  to  other  forms  including 
body  waves  that  can  propagate  into  the  half-space.  When  the  wavefield  reaches  the  second 
layered  structure,  modes  not  trapped  within  the  new  crustal  layer  will  leak  out  of  the  layer, 
rapidly  at  first,  then  at  a  steadily  decreasing  rate.  This  implies  that  wide  angle  reflections 
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of  energy  escaping  from  the  crustal  layer  towards  the  bottom  boundary  of  the  grid  could 
possibly  contaminate  transition  calculation  results.  However,  the  model  grids  do  not  extend 
far  enough,  in  the  x  direction,  beyond  the  transition,  for  this  to  be  a  problem.  The  energy 
will  encounter  the  rightmost  end  of  the  grid,  either  on  the  downgoing  or  the  upgoing  portion 
of  its  path,  rather  than  reaching  a  receiver  at  or  near  the  surface  as  a  wide  angle  reflection. 
Since  a  wavefront  which  has  a  small  angle  of  incidence  with  the  bottom  boundary  has  a  large 
angle  of  incidence  with  respect  to  the  end  boundary,  most  of  the  amplitude  of  the  wavefront 
incident  on  the  rightmost  end  boundary,  will  be  transmitted  rather  than  reflected.  Therefore, 
it  is  removed  from  the  grid.  Careful  grid  design  will  prevent  significant  contamination  from 
wide  angle  bottom  reflections. 

Changes  to  Lg  Wavetrains  on  Passage  through  a  C/O  Transition 
and  Their  Dependence  on  L 

When  a  wavefield  consisting  of  SH  L:  mode  sum  energy  passes  through  a  C/O  transition 
region  such  as  that  illustrated  in  figure  4a)  several  important  things  happen.  In  order  to 
clearly  explain  these  effects  and  their  variation  with  L  one  must  think  of  the  wavefield  in 
terms  of  several  components.  The  wavefield  may  be  divided  into  an  incident  wavefield.  a 
reflected  wavefield,  and  a  transmitted  wavefield.  The  incident  wavefield  is  composed  of  the 
trapped  modal  SH  Lg  energy  arriving  from  the  source,  and  is  present  and  complete  only  in 
the  region  of  continental  structure.  In  the  transition  region  some  components  of  the  incident 
wavefield  are  present,  but  as  the  wavefield  travels  through  the  transition  region  more  of  these 
components  are  converted  to  reflected  and  transmitted  wavefield  components.  The  reflected 
wavefield  is  composed  of  trapped  and  leaking  modes  which  have  been  reflected  from  the 
boundaries  in  the  transition  region  when  the  incident  wavefield  interacts  with  them.  The 
reflected  wavefield  travels  sourceward  through  the  crustal  layer,  of  the  transition  region, 
and  of  the  region  of  continental  structure.  The  transmitted  wavefield  is  composed  of  all 
the  energy  transmitted  through  the  transition  region  into  the  crustal  layer  of  the  region  of 
oceanic  structure,  and  all  the  energy  transmitted  into  the  mantle  layer.  The  transmitted 
wavefield  is  divided  into  four  components,  two  in  the  mantle  layer,  two  in  the  crustal  layer. 
The  forward  transmitted  wavefield  is  the  portion  of  the  incident  energy  which  crosses  the 
crust  mantle  boundary  (CMB)  into  the  mantle  layer  when  the  incident  wavefield  interacts 
with  the  CMB.  The  reverse  transmitted  wavefield  is  the  portion  of  the  reflected  energy 
which  crosses  the  CMB  into  the  mantle  layer  when  the  reflected  wavefield  interacts  with  the 
CMB.  The  forward  transmitted  wavefield  is  present  in  the  transition  region  and  the  region 


14 


of  oceanic  structure.  The  reverse  transmitted  wavefield  is  present  in  the  transition  region 
and  the  region  of  continental  structure.  The  directly  transmitted  wavefield  is  the  portion  of 
the  energy  in  the  incident  wavefield  that  passes  through  the  crustal  layer  of  the  transition 
region  and  into  the  crustal  layer  of  the  region  of  oceanic  structure  without  modal  conversion. 
The  transformed  transmitted  wavefield  is  the  component  of  the  incident  wavefield  which  is 
converted  into  a  new  set  of  trapped  and  leaking  modes  as  it  na«es  through  the  crustal  layer 
of  the  transition  region  into  the  crustal  layer  or  the  region  of  oceanic  structure. 

Each  of  the  components  of  the  wavefield  discussed  above  will  be  illustrated  in  the  time 
slices  or  seismograms  for  the  25  km  C/0  transition  region.  These  time  slices  and  seismo¬ 
grams  are  presented  as  an  example  of  the  general  transmission  properties  of  simple  transition 
regions.  The  variation  of  such  transmission  properties  with  L  will  also  be  illustrated.  The 
observed  properties  and  their  variations  with  L  will  then  be  explained. 

A  series  of  time  slices  for  the  25  km  C/O  transition  is  shown  in  Figure  6.  In  the  first 
time  slice  the  displacements  entering  the  grid  as  dark,  almost  vertical  bars  are  those  which 
produce  the  largest  peak  to  peak  amplitudes  in  the  seismograms.  The  third  time  slice  shows 
the  same  high  frequency  arrivals  after  they  have  propagated  about  halfway  through  the 
transition  region.  The  maximum  amplitudes  seen  within  the  transition  region  of  the  third 
time  slice  are  almost  twice  those  seen  in  the  unperturbed  layer  over  a  half-space  wavefield 
seen  in  the  first  two  time  slices.  So,  the  plotted  amplitudes  of  all  the  displacements  in  the 
third  time  slice  have  been  significantly  reduced  by  seeding,  making  the  same  disturbances 
appear  smaller.  In  the  fourth  time  slice  the  highest  amplitude  regions  of  the  wavefield 
have  propagated  past  the  rightmost  edge  of  the  grid.  The  maximum  amplitude  is  much 
smaller  so  the  same  disturbances  appear  to  have  much  larger  amplitudes.  The  fifth  time 
slice  shows  a  further  amplification  of  the  amplitude  of  the  disturbances  due  to  scaling.  The 
amplification  helps  make  the  highest  amplitude  portion  of  the  reflected  wavefield  (<  59f  of 
incident  amplitude)  visible  as  a  series  of  broken  vertical  bars  near  the  leftmost  end  of  the 
grid.  All  sections  of  the  incident  wavefield  show  triangular  regions  of  maximum  amplitude. 
The  extent  of  each  triangular  region,  in  the  x  direction,  increases  for  regions  of  the  wavefield 
incident  on  the  left  end  of  the  grid  at  a  later  time.  The  end  of  the  incident  wavefield  is  seen,  in 
the  fifth  time  slice,  about  midway  between  the  left  grid  edge  and  the  arrow,  to  the  right  of  the 
visible  reflected  component.  The  C/0  transmitted  wavefield  is  visible  in  the  last  three  time 
slices  as  disturbances  travelling  through  the  half-space  from  the  CMB  towards  the  rightmost 
edge  of  the  grid.  The  reverse  transmitted  wavefield  is  clearly  visible  in  the  fifth  time  slice. 
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and  present  in  the  last  three  time  slices,  as  disturbances  travelling  through  the  half-space 
from  the  CMB  towards  the  grid  bottom  and  the  source.  For  the  earliest  disturbances  in  the 
half-space  the  forward  transmitted  wavefield  propagates  along  a  path  approximately  parallel 
to  the  crust  mantle  interface  in  the  oceanic  region,  and  the  reverse  transmitted  wavefielc 
propagates  along  paths  nearly  parallel  to  the  CMB.  As  time  progresses  the  propagation 
paths  of  both  the  forward  and  reverse  transmitted  wavefields  approach  the  normal  to  the 
CMB.  Energy  escaping  from  the  transformed  transmitted  wavefield  as  it  propagates  through 
the  oceanic  crustal  layer  forms  a  forward  bending  arc  of  higher  amplitudes  in  the  mantle 
half-space.  As  the  hieh  amplitude  regions  of  the  forward  transmitted  wavefield  propagate 
longer  distances  or  propagate  along  paths  with  larger  vertical  components  the  tails  forming 
due  to  leakage  from  the  oceanic  crustal  layer  detach  from  the  from  the  forward  transmitted 
wavefield. 

Figure  7  shows  seismograms  recorded  at  the  receivers  along  a  surface  section.  These 
seismograms  illustrate  the  magnitude  of  the  concentration  of  amplitude  at  the  surface  of 
the  crust  seen  in  the  time  slices,  and  the  change  in  the  waveforms  as  the  wavefield  passer- 
through  the  transition.  Figure  8a  shows  the  variation  of  50  second  RMS  amplitude  with 
distance,  X2,  from  the  source  for  each  of  the  C/O  transition  FE  calculations.  Since  the 
relative  amplitudes  produced  using  the  different  C/O  transition  lengths  are  being  comparea 
the  absolute  values  of  the  amplitudes,  which  are  dependent  on  source  magnitude,  are  dot 
shown.  The  amplitudes  recorded  at  nodes  along  the  surface  of  the  crustal  layer  within  each 
transition  region  show  a  general  increase  as  one  moves  from  B  towards  D.  The  length  cf 
the  coda  with  amplitudes  above  one  third  of  the  maximum  peak  to  peak  amplitude  is  also 
increasing.  The  size  of  the  increase  in  RMS  amplitude  is  smallest  for  the  step  transition  and 
increases  as  L  increases.  Superimposed  on  the  general  increase  is  an  oscillatory  term.  The 
method  of  scaling  the  distance  coordinate  within  the  transition  region,  in  Figure  8.  makes 
the  coincidence  of  maxima  and  minima  at  approximately  the  same  fraction  of  the  transition 
length  for  all  L’s  considered  very  clear.  Increases  in  amplitude  before  reaching  distance  B 
are  due  to  energy  reflected  from  the  transition  boundaries  back  towards  the  source.  The 
amplitude  increase  near  B  is  largest  for  a  step  transition  and  decreases  as  L  increases.  The 
fluctuations  in  amplitude  following  the  transition  region  show  that  the  wavefield  has  no: 
stabilized  in  the  few  kilometers  beyond  the  transition  regions  illustrated  in  this  figure. 

Figure  9  shows  seismograms  recorded  on  depth  sections  at  distances  B  and  D.  Figure 
10  shows  the  variation  of  50  second  RMS  amplitude  with  depth  at  several  distances.  The 
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seismograms  at  the  depth  of  the  surface  of  the  oceanic  crustal  layer  show  a  5091  increase 
in  amplitude  as  the  energy  travels  from  B  to  D.  They  also  show  that  the  amplitude  of 
the  coda  becomes  a  larger  portion  of  the  maximum  peak  to  peak  amplitude.  The  nodes 
at  the  depth  of  the  bottom  of  the  oceanic  crustal  layer  show  a  small  decrease  in  RMS 
amplitude  across  the  transition.  Examination  of  the  nodes  at  B  and  D  with  depths  between 
these  pairs  (not  illustrated)  shows  that  the  increase  in  amplitude  is  largest  at  the  surface 
of  the  crust  and  decreases  rapidly  toward  the  base  of  the  crust.  The  amplitudes  of  the 
seismograms  transmitted  across  the  CMB  are  decreased  by  passage  through  the  transition 
region.  The  transmitted  waveforms  are  similar  to  the  incident  waveforms  and  show  an 
increasing  reduction  in  transmitted  amplitude  as  depth  increases.  Amplitude  is  seen  at 
depths  below  the  depth  of  the  bottom  of  the  continental  crustal  layer.  The  amplitudes  in 
the  mantle  layer  within  the  transition  region  are  maximum  immediately  below  the  CMB. 
and  decay  nearly  exponentially  with  depth  below  that  boundary.  The  amplitude  of  the 
transmitted  seismograms  recorded  at  nodes  below  the  oceanic  crustal  layer  and  above  the 
base  of  the  continental  crustal  layer,  at  D,  decreases  as  L  increases.  Below  this  depth  the 
trend  reverses.  At  distances  E  and  F,  further  from  the  transition  region,  the  magnitudes 
of  these  trends  decrease.  These  observations  are  further  support  that  as  L  increases  larger 
energies  are  transmitted  across  the  CMB  and  that  these  energies  propagate  downward  more 
rapidly. 

The  triangular  pattern  of  maximum  amplitudes  in  the  wavefields  (figure  6)  can  easily 
be  explained.  Consider  the  SH  Lt  energy  in  the  crustal  layer  a  superposition  of  of  post- 
critically  reflected  multiple  SH  wave  reflections.  The  critical  angle  is  about  51®.  Since  the 
wavefront  is  perpendicular  to  the  ray,  the  wavefronts  that  are  visible  as  the  triangular  regions 
of  maximum  amplitudes  should  show  the  observed  angles  of  incidence  with  the  boundary 
of  between  0®  and  39®.  The  increase  in  width  of  the  triangular  maxima  corresponds  to  an 
decrease  in  group  velocity  which  can  be  translated  to  a  decrease  in  period  and/or  a  larger 
contribution  from  higher  modes.  The  later  parts  of  the  seismogram  are  predominantly  higher 
mode  energy. 

Figure  1 1  illustrates  two  pairs  of  ray  diagrams,  one  for  a  25  km  C/O  transition  region,  the 
other  for  a  100  km  C/O  transition.  These  diagrams  are  useful  when  explaining  the  trends 
observed  in  the  results  presented  above,  and  their  dependence  on  L.  To  clearly  explain 
some  trends  multiple  reflections  within  the  transition  region  need  to  be  considered.  In  this 
paragraph  the  notation  and  relations  needed  in  such  explanations  will  be  defined.  For  the 
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case  of  multiple  reflections,  the  labeled  incident  angle  j,  will  become  j x,  and  the  inciden; 
angle  for  subsequent  reflections  within  the  transition  will  be  denoted  j„,  where  n  denotes 
the  n‘h  reflection  of  the  ray  after  it  enters  the  transition  region.  Similarly,  the  transmitted 
rays  for  the  ntk  reflection  of  a  ray  after  entering  the  transition  region,  will  be  denoted  j\. 
The  upper  transition  in  each  pair  shows  ray  paths  when  the  rays  encounter  the  CMB  of 
the  transition  before  the  COB.  When  these  rays  reach  the  CMB  their  incident  angles,  i. 
are  reduced  by  the  angle,  <pcm,  so  that,  j:  =  i  -  ocm.  As  long  as  the  ray  remains  in  the 
transition  region,  j„  =  i  -  n(<}>cm  +  ec«)  +  4>co-  The  lower  transition  in  each  pair  shows  the 
ray  paths  when  the  rays  encounter  the  COB  before  they  encounter  the  CMB.  For  this  case 
4>oe,  and,  <t>cm  are  interchanged  in  the  expression  for  j„.  The  energy  transmitted  across  the 
CMB,  with  jn  >  0,  forms  the  forward  transmitted  wavefield.  If  the  value  of  jn  becomes 
negative  before  it  reaches  the  end  of  the  transition  then  the  ray  has  been  reflected,  and 
successive  values  of  jn  will  become  more  negative  until  the  reflected  wavefront  reaches  the 
continental  structure.  The  energy  transmitted  across  the  CMB,  with  jn  <  0,  forms  the  reverse 
transmitted  wavefield.  The  angle  of  incidence  is  jae,  in  the  oceanic  crustal  layer,  and  jcm  for 
the  reflected  wave  in  the  continental  crustal  layer.  When  the  first  reflection  of  the  ray  exiting 
the  transition  region,  from  the  upper  or  lower  boundary  of  either  the  continental  or  oceanic 
crustal  layer,  is  the  nth  reflection  then  jctm  -  j„c  =  j„  -t-  4>»c,  or  jc<m  =  =  jn  +  <t>cm,  respectively. 

Also  useful  in  the  following  discussions  is  the  transmission  coefficient  as  a  function  of  the 
incident  angle,  j„.  For  the  examples  discussed  here  the  displacement  transmission  coefficient. 
T„(j„)  =  co50.'h°*6;Vo.(7'  )'  var>es  monotonically  between  zero  for  the  critical  angle  j„  =  51°,  and 
.76  for  normal  incidence,  j„  =  0. 

The  growth  in  amplitude  for  seismograms  at  crustal  surface  nodes  as  distances  range 
from  B  to  D  cat  be  explained  in  terms  of  energy  concentrated  in  the  thinning  crustal  layer 
of  the  transition  region.  As  j„  decreases  toward  zero  for  successive  n’s,  the  reflection  points 
at  the  boundaries  in  a  transition  region  are  separated  by  smaller  horizontal  distances.  This 
implies  that  the  density  of  rays  will  increase  as  the  crustal  layer  thins,  thereby  producing 
higher  amplitudes.  An  equal  concentration  would  be  expected  at  the  CMB  if  no  energy 
was  transmitted  across  that  boundary.  In  fact,  for  a  model  with  the  same  geometry  but 
with  the  mantle  half-space  replaced  by  water  to  disallow  transmission  across  the  CMB.  an 
amplification  of  75%  is  observed  at  both  boundaries.  However  for  continent  ocean  boundary 
models,  as  jn  and  thus,  j'n,  decrease  the  transmission  coefficient,  7V, ,  increases.  When  more 
energy  is  concentrated  at  the  CMB  more  energy  is  transmitted  across  it.  The  escape  of  enei g> 
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from  the  crustal  layer  to  form  the  transmitted  wavefield  also  explains  the  distribution  of 
amplitude  with  depth  within  the  crustal  layer  of  the  transition.  The  amount  of  concentration 
of  amplitude  increases  as  the  distance  from  the  CMB  increases  and  the  effects  of  the  energy 
escaping  across  it  become  weaker. 

Extending  the  above  arguments  explains  the  increase  in  the  amount  of  concentrated 
energy,  the  increase  in  the  magnitude  of  the  transmitted  wavefield.  and  the  decrease  in 
the  amplitude  of  the  reflected  wavefield  as  L  increases.  As  L  increases  the  angles  o and 
C'o  decrease  causing  j„  to  decrease  more  slowly  and  allowing  the  number  of  reflections,  m. 
within  a  transition  to  increase.  When  L  increases,  the  increase  in  m  causes  an  increase  in  ray 
density  and  explains  the  increase  in  concentrated  amplitude  at  the  COB.  When  L  increases 
the  amount  of  energy  concentrated  at  the  CMB  also  increases.  The  transmission  of  this 
energy  across  the  CMB  causes  the  amplitude  of  the  transmitted  wavefield  to  increase  as  L 
increases.  Since  the  amount  of  amplitude  concentrated  at  the  COB  and  thus  transmitted 
into  the  oceanic  crustal  layer,  and  the  amount  of  energy  concentrated  at  the  CMB  and  thus 
transmitted  across  the  CMB  both  increase  as  L  increases  the  amplitude  of  the  reflected 
wavefield  must  decreases  as  L  increases.  Let  seismogram  A  be  a  hybrid  seismogram  recorded 
on  depth  section  A  (Figure  4)  in  a  C/O  transition  calculation.  Let  seismogram  B  be  the 
hybrid  seismogram,  recorded  during  the  continental  reference  calculation,  at  the  same  depth 
and  distance  from  the  source  as  seismogram  A.  The  reflected  component  of  is  measured  by 
analyzing  the  component  which  remains  when  seismogram  B  is  subtracted  from  seismogram 
A.  Reflected  amplitudes  measured  in  this  manner  show  a  clear  decrease  as  L  increases. 

Figure  1 1  can  be  used  to  explain  how  and  why  the  direction  of  propagation  of  the  forward 
and  reverse  transmitted  wavefields  varies  as  time  progresses.  When  the  arrivals  of  the  highest 
group  velocity  reach  the  CMB  the  resulting  forward  transmitted  wavefield  travels  along  a 
path  approximately  parallel  to  the  crust  half-space  interface  in  the  oceanic  structure,  and 
the  resulting  reverse  transmitted  wavefield  travels  along  a  path  almost  parallel  to  CMB 
The  incident  angles  j„  are  maximum.  As  |y„|  for  a  given  n,  and  the  group  velocity,  decrease 
the  values  of  j'„  decrease.  Thus,  the  propagation  direction  of  both  transmitted  wavefields 
approach  the  normal  to  the  CMB.  As  L  increases  larger  values  of  n  are  possible  within 
the  transition  region.  For  each  value  of  n  a  transmitted  wavefield  is  created  at  the  CMB. 
Therefore,  the  range  of  propagation  directions  seen  at  a  given  time  increases  as  L  increases. 

The  directly  transmitted  wavefield  and  the  transformed  transmitted  wavefield  can  also  be 
easily  explained.  For  reflections  points  immediately  preceding  the  beginning  of  the  transition 
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region  sonic  rays  (*  >  90  —  arctan{27 j L))  can  propagate  through  the  transition  region  without 
interacting  with  either  the  CMB  or  the  COB.  These  rays,  which  form  the  directly  transmitted 
wavefield.  continue  to  bounce  with  the  same  post-critical  angles  of  incidence  at  the  crust 
mantle  interface  and  the  free  surface  as  they  did  in  the  continental  structure.  They  represent 
the  energy  transmitted  through  the  transition  directly,  without  changes  in  period  or  phase 
velocity.  For  these  rays  no  mode  conversion  has  occurred  and  no  energy  has  been  converted 
to  other  phases.  These  rays  form  the  unchanged  components  of  the  eigenfunctions  defining 
each  mode.  The  range  of  angles,  i,  over  which  the  direct  transmission  occurs  decreases  as 
L  increases.  Thus,  more  unconverted  energy  is  transmitted  through  shorter  transitions.  In 
contrast,  the  transformed  transmitted  wavefield  is  composed  of  rays  which  interact  with 
the  CMB  or  COB  before  being  transmitted  into  the  oceanic  crustal  layer.  Such  rays  show 
changes  in  group  velocity  that  indicate  modal  conversions  are  occurring.  Energy  is  also 
converted  into  other  phases  as  it  escapes  across  the  CMB.  As  discussed  above,  the  energies 
transmitted  into  the  oceanic  crustal  layer  and  across  the  CMB  both  increase  as  L  increases. 
As  L  increases,  this  increased  transmitted  energy  is  partially  counteracted  by  the  effects  of 
decreased  transmission  of  unconverted  energy.  This  helps  explain  why  the  differences  due 
to  L  are  small.  It  should  be  noted  that  increasing  L  reduces  the  rate  at  which  the  crustal 
thickness  changes,  so  at  a  given  distance  from  the  beginning  of  a  transition  region  the  change 
in  amplitude  transmitted  across  the  CMB  per  unit  horizontal  distance  becomes  smaller  as 
L  increases.  However,  integrated  over  the  entire  length  of  the  transition  region  the  amount 
of  energy  transmitted  across  the  CMB  increases. 

The  diagrams  in  Figure  11  demonstrate  why  energy  crosses  the  crust  mantle  interface 
in  the  oceanic  structure  following  the  transition,  why  this  transmitted  energy  is  maximum 
near  the  transition  and  decreases  as  the  distance  from  the  transition  increases,  and  why 
th  is  transmitted  energy  increases  as  L  decreases.  The  paths  illustrated  in  the  second  25 
km  transition  outline  show  one  way  modes  which  are  not  of  appropriate  frequencies  to 
be  trapped  within  the  oceanic  crustal  layer,  enter  that  layer.  Analogous  paths  for  rays 
which  interact  first  with  the  CMB  or  rays  which  interact  multiple  times  with  the  crustal 
boundaries  in  the  transition  region  also  exist.  All  such  rays  produce  rays  at  pre-critical. 
joc  <  51°.  angles  within  the  oceanic  crust.  For  any  such  ray,  the  transmission  coefficient  at 
each  successive  reflection  from  the  crust  mantle  interface  is  T(joe)-  Thus,  the  proportion  of 
reflected  amplitude  remaining  in  the  crust  after  the  n"1  bounce  is  ,4„  =  (1  -  T(j0 c))°.4.  where 
A  is  the  original  amplitude.  Clearly,  the  amount  of  escaping  energy  decreases  with  distance. 
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I-'oi  a  shorter  transition  jn,  for  a  given  i,  is  smaller.  Thus,  a  larger  range  of  angles  joc  <  51° 
are  produced  for  each  path  through  the  shorter  transition.  This  means  a  higher  density 
of  pre-critical  rays  and  higher  amplitudes  in  and  escaping  from  the  oceanic  crustal  layer  of 

shot  ter  transitions. 

The  use  of  the  interpretation  of  Ls  waves  in  terms  of  multiply  reflected  SH  rays  to  the 
interpret  the  FE  results  has  been  instructive.  However,  the  many  interacting  effects  involved 
in  determining  the  properties  of  the  reflected,  converted,  and  transmitted  wavefieids  would 
be  extremely  difficult  to  predict  using  this  approach.  Thus,  the  FE  method  is  necessary  to 
determine  which  effects  are  important  and  which  need  not  be  considered. 

Changes  to  L}  Wavetrains  on  Passage  through  a  O/C  TVansit  ion  Region 

and  Their  Dependence  on  L 

The  passage  of  a  wavefield  consisting  of  SH  type  L,  mode  sum  energy  for  a  continental 
laver  over  a  half-spac  model,  which  has  passed  through  a  C/O  transition  of  50  km  length, 
through  a  O/C  transition  such  as  that  illustrated  in  Figure  4b'>  has  several  effects  on  the 
wavefield.  Results  of  the  calculations  for  propagation  through  a  50  km  O/C  transition  will 
be  presented  in  detail  as  to  illustrate  these  effects.  The  O/C  transition  tests  using  a  variety 
of  L's  will  be  summarized  to  illustrate  how  these  effects  depend  on  L.  The  noted  effects  and 
their  variations  with  L  will  be  explained.  Before  proceeding  to  these  discussions  the  errors 
introduced  due  to  the  truncation  inherent  in  coupling  the  C/O  transition  FE  results  into  the 
O/C  transition  grid  should  be  mentioned.  The  finite  vertical  extent  of  nodes  driven  by  forcing 
functions  at  the  leftmost  edge  of  the  O/C  transition  grid  will  cause  a  vertical  truncation  of 
the  incoming  wavefield.  Also,  any  reflections  included  in  the  seismograms  recorded  to  use  as 
forcing  functions  will  be  added  to  the  forward  propagating  wavefield  in  the  O/C  calculation. 
These  effects  tend  to  increase  the  amplitudes  at  the  surface  of  the  oceanic  layer  near  the  left 
end  of  the  grid.  The  uncertainties  introduced  by  the  coupling  process  appear  to  increase  the 
amplitudes  in  the  second  grid  by  as  much  as  3%.  Therefore,  they  could  possibly  lead  to  a 
slight  underestimate  of  the  magnitude  of  the  attenuation  effect,  but  should  not  lead  to  an 
overestimate. 

A  series  of  time  slices,  shown  in  Figure  12,  illustrate  the  FE  calculation  results  for  the 
50  km  long  O/C  transition.  The  highest  amplitude  concentrations  visible  in  first  four  slices 
correspond  to  the  maximum  amplitudes  in  the  seismograms.  The  first  two  time  slices  illus¬ 
trate  the  wavefield  travelling  through  the  oceanic  structure.  The  third  time  slice  illustrates 
the  passage  of  the  highest  amplitude  portions  of  the  wavefield  through  the  transition  region 
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In  the  fourth  time  slice  these  highest  amplitude  regions  appear  at  the  rightmost  end  of  the 
grid.  In  the  fifth  time  slice  they  have  exited  the  grid.  The  maximum  amplitude  in  a  given 
time  slice  decreases  with  time  in  the  third  to  fifth  time  slice.  The  maximum  amplitudes  in 
the  transition  region  decrease  as  the  wavefront  proceeds  through  the  transition  region.  Thus, 
normalizing  by  the  largest  amplitude  in  each  slice  causes  identical  disturbances  to  increase 
in  extent  and  intensity  in  successive  slices.  Triangular  regions  of  maxima  can  be  seen  in 
the  oceanic  crust  but  they  are  not  nearly  so  clear  as  those  seen  in  the  continental  crust  of 
the  C/O  transition.  The  third  and  fourth  time  slices  clearly  show  that  the  energy  which 
was  previously  trapped  in  the  oceanic  crustal  layer  distributes  itself  throughout  the  crus¬ 
tal  layer  as  it  passes  through  the  transition  region  and  the  region  of  continental  structure. 
The  amplitudes  remain  much  larger  near  the  surface,  and  decay  rapidly  with  depth.  The 
fifth  time  slice  shows  triangular  regions  maximum  amplitude  are  beginning  to  appear  in  the 
continental  structure.  Additional  calculations  show  further  propagation  in  x  is  necessary  to 
stabilize  the  wavefield  in  the  region  of  continental  structure  . 

The  disturbances  seen  in  the  half-space  of  the  oceanic  region  of  the  0/C  transition  grid 
contain  significant  energy.  Examining  the  time  slices  shows  that  for  a  subset  of  these  distur¬ 
bances,  equivalent  to  the  forward  transmitted  wavefield  discussed  earlier,  the  z  component 
of  translation  in  the  half-space  is  increasing  for  successive  groups  of  disturbances.  This  is 
due  to  the  interactions  at  the  CMB  in  the  C/O  transition  as  explained  earlier.  It  is  more 
visible  in  these  time  slices  than  the  C/O  transition  time  slices  due  to  the  increased  distance 
from  the  CMB.  It  is  clear  that  these  packets  of  amplitude  are  translated  rapidly  enough  in  z 
that  some  will  pass  into  the  half-space  below  the  continental  crustal  layer  before  they  reach 
the  O/C  transition  region.  Thus,  the  energy  contained  within  them  escapes  the  system  and 
is  not  reconverted  to  L,  energy  when  the  wavefield  passes  through  the  O/C  transition.  The 
amount  of  energy  escaping  from  the  system  in  this  manner  increases  as  L,  or  the  length  of 
the  intermediate  oceanic  path,  increases. 

Figure  13  shows  a  surface  section  of  seismograms  for  a  50  km  O/C  transition.  Figure  8b 
shows  the  variation  of  50  second  RMS  amplitude  with  distance,  X2,  from  the  source  for  each 
of  the  O/C  transition  FE  calculations.  A  decrease  in  RMS  amplitudes  is  seen  as  one  moves 
along  the  surface  of  the  oceanic  crustal  layer  approaching  the  transition  region.  This  decrease 
is  clear  for  short  RMS  windows  and  almost  vanishes  for  window  length  of  65  sec  or  more. 
This  indicates  that  the  portion  of  the  coda  with  amplitudes  comparable  lo  the  maximum 
amplitude  is  decreasing  in  length.  Superimposed  on  this  slow  decrease  in  amplitude  is  an 


22 


oscillatory  term  with  amplitude  approximately  15%  of  the  mean  amplitude  at  the  surface 
in  the  oceanic  structure.  The  RMS  amplitudes  recorded  at  nodes  along  the  surface  of  the 
crustal  layer  within  each  0/C  transition  region  decrease  as  one  moves  from  B  towards  D, 
The  size  of  this  decrease  is  smallest  for  the  step  transition  and  increases  as  L  increases.  The 
rate  of  decay  shows  some  tendency  to  become  smaller  towards  the  end  of  the  transition.  The 
superimposed  oscillations  seen  in  the  O/C  transition  results  are  much  smaller  than  those 
seen  in  the  C/O  transition  results.  They  appear  to  be  confined  near  the  beginning  of  the 
transition  region.  The  oscillation  in  amplitude  following  the  transition  region  shows  that  the 
wavefield  has  not  stabilized  in  the  few  kilometers  beyond  the  transition  regions  illustrated 
in  figure  Sb.  Additional  results  which  are  not  illustrated  show  the  RMS  amplitude  continues 
to  decrease  slowly  until  the  wavefield  readjusts  to  the  new  continental  structure. 

Figure  14  shows  the  variation  of  50  second  RMS  amplitude  with  depth.  At  the  depth 
of  the  surface  of  the  oceanic  crustal  layer  propagation  through  the  transition  reduces  the 
RMS  amplitude  by  more  than  a  factor  of  two,  and  the  peak  to  peak  amplitude  by  25%.  This 
indicates  that  the  length  and  amplitude  of  the  coda  are  decreasing  more  rapidly  than  the 
amplitude  of  the  largest  arrival.  At  depths  between  the  surface  of  the  continental  crust  and 
the  surface  of  the  oceanic  crust  amplitude  decreases  with  depth.  Then  energy  incident  on 
the  MCB  is  largely  transmitted  across  the  MCB  into  the  crustal  layer.  Some  of  this  energy 
eventually  escapes  across  the  CMB  in  the  region  of  continental  structure.  The  amplitudes 
below  the  crustal  layer  of  the  transition  region  decrease  slowly  with  depth  or  even  increase 
with  depth.  An  increase  of  amplitude  with  depth  below  the  crustal  layer  indicates  that 
energy  is  travelling  down  toward  the  grid  bottom  and  out  of  the  system  being  considered. 

Most  of  the  energy  incident  on  the  MCB  is  transmitted  across  it.  Thus,  in  the  last  four 
time  slices,  the  reflections  from  the  MCB  are  not  easily  separated  from  the  later  portions 
of  the  incident  wavefield.  The  pattern  of  displacements  seen  in  the  oceanic  half-space  is 
distorted  when  the  energy  producing  it  is  transmitted  across  the  MCB.  The  uppermost 
portion  of  a  disturbance  crosses  the  boundary  first  and  begins  to  move  more  slowly.  The 
remainder  of  the  disturbance,  still  in  the  mantle  layer  continues  to  move  with  the  mantle's 
faster  velocity.  As  the  disturbance  moves  through  the  length  of  the  transition,  an  increasing 
section  of  it  moves  into  the  crustal  layer.  This  results  in  a  slope  being  superimposed  on 
the  portion  of  the  disturbance  which  has  propagated  back  into  the  crustal  layer.  The  slope. 
Slope  -  is  dependent  on  L,  the  velocity  in  the  crustal  layer,  vc,  and  the  velocity  in 

the  mantle,  \>M .  This  slope,  in  this  case  about  tan  52°,  is  observed  easily  in  both  the  second 
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and  third  time  slices.  It  is  also  seen  in  the  fourth  and  fifth  time  slices.  However,  in  these 
time  slices  the  amplitude  incident  from  the  oceanic  crustal  layer  dominates  that  transmitted 
from  the  oceanic  mantle  layer  making  observation  of  the  slope  more  difficult. 

The  effect  of  the  O/C  transition  region  on  the  energy  previously  travelling  in  the  oceanic 
crustal  layer  will  now  be  discussed.  Concentrations  of  amplitude  incident  from  the  oceanic 
crustal  layer  travel  up  the  surface  of  the  OCB.  As  they  propagate  up  the  OCB  the  lower 
edge  of  the  amplitude  concentrations  are  no  longer  constrained  by  the  lower  edge  of  the 
oceanic  crustal  layer.  Consequently  the  energy  can  migrate  downwards  towards  the  depth  of 
the  base  of  the  continental  crustal  layer.  This  diffusion  of  energy  can  be  explained  using  a 
mechanism  which  is  the  converse  of  that  used  to  explain  the  concentration  of  energy  in  the 
C/O  transition.  The  incident  angle  for  the  ray  perpendicular  to  each  wavefront.  j„,  increases 
as  n  increases.  Subsequent  reflection  points  are  separated  by  larger  horizontal  distances 
producing  a  dilution  in  ray  density  and  thus  a  decrease  in  amplitude.  However,  there  is  an 
important  difference  that  simplifies  the  analysis  of  the  O/C  transition.  Incident  angles.  jx. 
for  all  rays  trapped  in  and  travelling  through  the  oceanic  crustal  layer  are  greater  than  51°. 
For  n  >  1  these  angles  increase.  Therefore,  no  energy  originally  trapped  in  the  oceanic  crust 
is  transmitted  across  the  MCB,  or  reflected  back  toward  the  source. 

Observations  indicate  that  the  length  of  the  intermediate  oceanic  path  between  the  C/O 
and  O/C  transitions  is  important.  The  results  of  this  study  suggest  some  reasons  why  this  is 
so.  For  short  paths  energy  transmitted  in  to  the  mantle  layer  has  little  time  to  travel  towards 
the  bottom  of  the  grid  and  thus  out  of  the  region  of  consideration  before  much  of  it  passes 
back  into  the  crustal  layer  in  the  O/C  transition  region.  This  implies  that  there  may  be  a 
critical  length  of  intermediate  oceanic  path  beyond  which  enough  of  the  energy  has  escaped 
from  the  region  of  interest  that  amplitudes  of  the  attenuated  L,  recorded  after  the  O/C 
transition  would  be  reduced  sufficiently  to  explain  the  observed  data.  To  completely  analyze 
this  assertion  requires  the  use  of  the  RT  coupling  method  for  continuing  FE  calculations 
through  a  plane  layered  structure  using  RT  integration  and  PM  Green’s  functions.  The 
accuracy  of  the  numerical  implementation  of  this  method  must  be  established  before  such 
calculations  can  be  presented.  Therefore,  the  discussion  of  the  numerical  implementation  of 
this  method  and  of  the  results  examining  the  effects  of  the  intermediate  path  length  in  the 
oceanic  structure  will  be  left  for  other  papers. 
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Conclusions 


The  effects  produced  when  a  SH  Lt  mode  sum  wavefield  is  propagated  through  a  C/O 
transition  region,  and  their  dependence  on  L,  will  now  be  summarized.  Amplitudes  at  the 
surface  of  and  in  the  crust  are  amplified  by  as  much  as  50%  as  the  wavefield  passes  through 
the  transition  region.  These  amplitude  increases  in  the  crustal  layer  have  been  explained  bv 
the  increased  ray  density  in  the  thinning  crustal  layer.  The  amount  of  amplification  increases 
as  L  increases.  Increased  numbers  of  reflections  within  a  longer  transition  region,  cause  an 
increase  in  ray  density,  explaining  the  increase  in  energy  as  L  increases.  Amplification  is 
maximum  at  the  surface  and  decreases  with  depth  until  the  base  of  the  oceanic  crust  is 
reached.  Amplitude  is  transmitted  across  the  CMB  and  but  not  across  the  COB.  Thus,  the 
amplification  is  maximum  at  the  COB  and  decays  as  the  CMB  is  approached  and  the  effects 
of  energy  transmitted  across  the  CMB  increase.  The  amount  of  energy  transmitted  across  the 
CMB  increases  as  L  increases.  The  translation  direction  of  the  forward  transmitted  wavefield 
has  a  larger  vertical  component  later  in  the  seismogram  when  it  approaches  the  normal  to 
the  CMB.  Thus,  for  a  longer  transition,  with  a  smaller  4>cm  and  thus  a  more  horizontal  CMB. 
the  path  length  in  the  oceanic  structure  will  have  a  larger  vertical  component.  This  will 
allow  the  transmitted  energy  to  escape  the  system  more  readily  as  L  increases. 

After  the  wavefield  has  travelled  through  the  transition  region  additional  energy  escapes 
from  the  crustal  layer  in  the  oceanic  structure.  This  leakage  is  largest  near  the  end  of 
the  transition  region  and  decreases  with  distance  away  from  it.  The  size  of  the  decrease 
in  amplitude  is  controlled  by  the  energy  in  wavefronts  propagating  through  the  oceanic 
crustal  layer  with  pre-critical  angles  of  incidence.  At  each  successive  reflection  of  such  a 
wavefront  at  the  crust  half-space  interface  in  the  oceanic  structure,  a  fraction  of  the  total 
amplitude  remaining,  T(joe)  is  transmitted  across  the  interface.  Thus,  the  total  amplitude 
of  this  component  decays  most  rapidly  in  the  first  kilometers  of  the  oceanic  structure  while 
the  incident  wave  still  has  significant  amplitudes. 

Next,  a  summary  of  the  important  observations  for  the  O/C  transition  will  be  given. 
As  the  wavefield  travels  through  the  O/C  transition  region  the  surface  amplitude  decreases 
rapidly.  The  energy  forming  the  concentrations  of  amplitude  previously  trapped  in  the 
oceanic  crustal  layer  travels  up  the  OCB  while  spreading  itself  throughout  the  crustal  layer. 
The  maxima  in  these  regions  remain  at  the  surface.  The  amplitude  decreases  rapidly  with 
depth  down  to  the  crust  mantle  interface.  The  amplitude  incident  upon  the  MCB  from 
the  oceanic  half-space  is  transmitted  across  it  with  a  resulting  distortion  in  the  wavefield. 
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The  form  of  this  distortion  is  the  superposition  of  a  slope,  Slope  =  fan-1  [ — * ~ 1  - ' ] .  onto  the 
disturbance  incident  at  the  MCB  after  it  has  passed  into  the  continental  crustal  layer. 

Finally,  attenuation  of  La  propagating  along  a  mixed  path  with  simple  transition  regions 
and  an  oceanic  path  length  of  147  km,  is  not  sufficient  to  explain  the  observed  values.  The 
highest  attenuations  observed  produced  a  reduction  of  a  factor  of  two  to  three  in  ampli¬ 
tude.  Observations  of  the  trends  due  to  varying  L  indicate  that  a  longer  C/O  transition 
could  increase  the  attenuation.  But,  even  using  the  optimal  L  attenuation  of  more  than 
a  factor  of  four  would  not  be  expected.  Thus,  other  factors  and/or  more  complex  struc¬ 
tures  must  be  considered  to  explain  the  observed  attenuation  of  Lg  along  mixed  paths. 
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TABLE  1 

GRID  CHARACTERISTICS  FOR  TRANSITION  MODELS 
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Figure  1:  Types  uf  models  used  in  studies  of  propagation  of  Love  waves  across  continent 
ocean  boundaries,  in  order  of  increasing  complexity,  a)  show's  two  homogeneous  layered 
regions.  1  and  3,  separated  by  an  intermediate  region,  2,  in  which  elastic  properties  van 
smoothy  between  their  values  in  regions  1  and  2.  b)  shows  a  layer  over  a  half-space  with 
a  step  change  in  the  thickness  of  the  layer,  c)  shows  a  model  with  a  smooth  change  in 
thickness,  either  at  the  surface  or  the  Moho.  d)  shows  a  model  with  a  smooth  change  m 
thickness  both  at  the  surface  and  the  Moho.  The  variations  in  a),  c),  and  d)  occur  in  a 
transition  region  of  length  L. 

Figure  2:  Geometry  used  to  explain  the  coupling  of  PM  seimograms  from  a  source 
outside  a  FE  grid  into  a  FE  grid.  The  two  long  horizontal  lines  show  the  free  surface  and 
the  boundary  between  the  layer  and  the  half-space.  The  source  is  shown  as  an  asterisk.  Two 
columns  of  FE  nodes  are  shown  as  dots.  The  vertical  line  connecting  the  dots  and  the  short 
horizontal  line  perpendicular  to  it  are  grid  edges.  The  heavy  solid  line  from  the  source  to  the 
surface  receiver  in  the  grid  denotes  the  direct  analytic  seismogram,  the  solid  lines  between 
the  source  and  the  grid  edge  nodes  denote  the  direct  forcing  functions,  and  the  dotted  lines 
indicate  the  source  to  receiver  paths,  for  the  sources  created  by  the  application  of  the  forcing 
functions  at  the  grid  edge.  Integration  over  all  dashed  paths  gives  the  hybrid  seismogram 

Figure  3:  Sample  PM  mode  sum  synthetic  seismogram  for  SH  Lg  1500  km  from  the 
source.  The  seismogram  includes  modes  with  periods  between  0.5  sec  and  100  sec.  The 
fundamental  and  the  first  five  higher  modes  are  used.  The  seismogram  has  been  band  pass 
filtered  between  0.01  and  1  hz,  and  has  the  WWSSN  short  period  instrument  has  been 
convolved  with  it.  The  box  labeled  input  window  shows  the  time  window  used  to  select  the 
forcing  function.  The  numbers  along  the  axis  indicate  the  group  velocities  of  the  arrivals. 

Figure  4:  Explanation  of  terms  and  illustration  of  the  model  classes  used  to  describe  the 
behavior  on  passage  through  a  transition  region.  The  heavy  line  betu-een  the  water  layer 
and  the  crustal  layer  is  the  surface.  The  sloping  portion  of  this  surface  is  the  continent  to 
^  ocean  boundary  (COB)  for  the  C/O  transition  model  and  the  ocean  to  continent  boundary 

(OCB)  for  the  O/C  transition  model.  Similarly,  the  sloping  dashed  line  between  the  crust 
and  mantle  layers  is  the  crust  to  mantle  boundary  (CMB1  for  the  C/O  transition  and  the 
mantle  to  crust  boundary  (MCB)  for  the  O/C  case.  The  length  of  transition.  L.  is  the 
distance  from  B  to  D,  B  is  referred  to  as  the  beginning  of  the  transition,  D  as  the  end  of  the 
transition,  C  as  the  center  of  the  transition.  A  is  5  km  from  B,  E  is  5  km  from  D. 

Figure  5:  Seismograms  recorded  at  the  surface  for  X2=1520  and  1540  km.  The  first 
seismogram  in  each  group  is  a  direct  synthetic,  the  second  is  a  hybrid  synthetic  determined 
using  the  long  grid  (200x100),  the  third  and  fourth  are  hybrid  synthetics  determined  using 
a  short  (100x100)  grid  with  and  without  the  BC  respectively.  The  number  beside  each  of 
the  hybrid  synthetics  is  the  RMS  amplitude  ratio  between  that  hybrid  seismogram  and  the 
direct  synthetic.  Seismograms  have  been  band  pass  filtered  between  .01  and  1  Hz  and  have 
been  convolved  with  the  short  period  WWSSN  instrument.  Arrival  times  of  the  single  and 
multiple  reflections  are  indicated  by  the  arrows  below  each  group  of  seismograms. 

Figure  6:  Time  slices  recorded  during  the  25  km  C/O  transition  calculation.  The  crustal 
layer  and  the  ends  and  bottom  of  the  FE  grid  which  contain  it  are  outlined.  A  larger 
displacement  produces  a  larger  point.  Thus,  high  displacement  areas  appear  darker.  Sixty 
forcing  functions  were  used.  The  time  since  the  initiation  of  the  FE  calculation  is  shown 
above  the  right  end  of  each  grid.  The  arrow  above  a  time  slice  shows  the  location  to  which 
the  disturbances  moving  at  3.5  km/sec,  seen  at  the  left  end  of  the  previous  time  slice,  have 
moved  in  the  intervening  duration. 
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Figure  7-  Surface  section  oi  seismograms  recorded  for  a  25  km  C/O  transition  calcuiat  ioi: 
Hie  numbers  to  the  left  of  each  seismogram  indicate  the  location  of  the  node  at  which  that 
seismogram  is  recorded.  Seismogram  1  was  recorded  at  B.  seismogram  11  at  D  The  numhn< 
increase  one  moves  farther  from  the  source,  1  he  spacing  between  receivers  in  the  transition 
70 ne  is  uniform.  The  numbers  above  the  right  end  of  each  seismogram  give  die  60s  RMS 
amplitude  of  that  seismogram. 

Figure  S:  Variation  in  RMS  surface  amplitude  with  distance.  X2.  from  the  source  for 
each  C/0  transition  calculation  (upper  sketch.  60  second  RMS),  and  for  each  0/C  transition 
calculation  (lower  sketch,  50  second  RMS).  RMS  amplitude  is  the  y  coordinate.  The  two 
verdcal  lines,  labeled  B  and  D,  indicate  the  limits  of  the  transition  region.  Outside  the 
transition  region  the  scale  is  uniform.  The  distance  scale  (x  coordinate)  in  the  transition 
region  is  different  for  each  L.  Within  the  transition  region  all  distances  are  plotted  with 
respect  to  an  origin  at  B,  as  fractions  of  L.  Each  line,  labelled  with  the  L  it  represents,  was 
generated  using  RMS  measurements  at  intervals  of  5  km  along  the  crustal  surface. 

Figure  0:  Depth  sections  of  seismograms  recorded  at  B  and  D  for  a  25  km  C/O  transi 
lion  calculation.  The  numbers  above  the  right  end  of  each  seismogram  show  the  6Us  RMS 
amplitude.  Successive  rows  show  pairs  of  receivers  at  increasing  depths.  The  first  row  of 
seismograms  shows  receivers  at  the  depth  of  the  surface  of  the  continental  crust,  the  second 
row  receivers  at  the  depth  of  the  surface  of  the  oceanic  crust,  the  third  row  receivers  at  the 
depth  of  the  base  of  the  oceanic  crust,  the  fifth  row  receivers  at  the  depth  of  the  base  of  the 
continental  crust,  and  the  fourth  row  receivers  midway  between  the  third  row  and  the  fifth 
row. 

Figure  10:  Variation  of  55  second  RMS  amplitude  with  depth  on  depth  sections  A.  B.  C’. 
D.  E.  of  Figure  4.  and  F,  located  10  km  beyond  E,  for  each  C/O  transition  model.  Each  plot 
illustrates  results  for  one  depth  section.  The  letter  at  the  upper  left  of  each  plot  identifies 
the  location  of  the  depth  section.  Amplitudes  for  the  surface  node  and  nodes  equally  spaced 
down  the  depth  section  (Az  =  2.5  km)  are  shown.  The  three  solid  horizontal  lines  are. 
from  top  to  bottom,  the  surface  and  base  of  the  oceanic  crustal  layer,  and  the  base  of  the 
continental  crustal  layer.  The  dashed  horizontal  line  shows  the  depth  of  the  base  of  the 
crustal  layer  in  depth  section  C,  at  the  center  of  the  transition  region. 

Figure  11:  Illustration  of  ray  paths  within  a  25  km  and  a  100  km  C/O  transition  struc¬ 
ture.  Diagrams  are  drawn  to  scale.  Each  diagram  shows  rays  with  angles  of  incidence,  i.  55°. 
65°,  and  75°,  at  the  free  surface  and  crust  mantle  interface  in  the  continental  portion  of  the 
model.  These  angle  sample  the  possible  range  of  post-critical  incident  angles  (>  5l/circ).  The 
incident  angle,  jn,  and  the  transmission  angle.  j'n ,  at  the  crustal  boundaries  in  the  transition 
regions  are  labeled.  The  upper  diagram  in  each  pair  shows  propagation  paths  for  rays  that 
encounter  the  CMB  before  the  COB.  The  lower  diagram  shows  propagation  paths  for  rays 
that  encounter  the  COB  before  the  CMB.  The  arrows  indicate  the  direction  of  propagation 
of  the  wavefront  along  the  ray. 

Figure  12:  Time  slices  recorded  during  the  50  km  O/C  transition  calculation.  Details 
are  identical  to  Figure  7.  40s  RMS  amplitudes  are  shown. 

Figure  13:  Surface  section  of  seismograms  recorded  for  a  50  km  O/C  transition  calcula¬ 
tion.  Details  explained  in  Figure  8. 

Figure  14:  Variation  of  55  second  RMS  amplitude  with  depth.  Details  are  as  explained 
in  Figure  10. 
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Modeling  regional  Love  Waves : 
imperial  Valley  to  Pasadena 

Phyllis  Ho-Liu  and  Donald  V.  Helmberger 

Abstract 

Long  period  seismograms  recorded  at  Pasadena  of  earthquakes  occurring  along  a  profile  to 
Imperial  Valley  are  studied  in  terms  of  source  phenomena  versus  path  effects.  Some  of  the  events 
have  known  source  parameters,  determined  by  teleseismic  or  near-field  studies,  and  are  used  as 
master  events  in  a  forward  modeling  exercise  to  derive  the  Green's  functions  (SB  displacements 
at  Pasadena  due  to  a  pure  strike-slip  or  dip-slip  mechanism)  that  describe  tha  propagation  effects 
along  the  profile.  Both  timing  and  waveforms  of  records  are  matched  by  synthetics  calculated 
from  2-dimensional  velocity  models.  The  best  2-dimensional  section  begins  at  Imperial  Valley 
with  a  thin  crust  containing  the  basin  structure  and  thickens  towards  Pasadena.  The  detailed 
nature  of  the  transition  zone  at  the  base  of  the  crust  controls  the  early  arriving  shorter  periods 
(strong  motions)  while  the  edge  of  the  basin  controls  the  scattered  longer  period  surface  waves. 
Shallow  events  in  the  basin  are  easily  distinguished  from  deep  events  as  well  as  the  amount 
of  strike-slip  versus  dip-slip  motions  Involved  from  the  waveform  characteristics  alone.  Those 
events  rupturing  the  sediments,  such  as  the  1979  Imperial  Valley  earthquake,  can  be  recognized 
easily  by  a  late  arriving  scattered  Love  wave  which  has  been  delayed  by  the  very  slow  path 
across  the  shallow  valley  structure. 

Introduction 

The  major  objectives  of  this  paper  are  to  initiate  a  field  of  study  which  will  allow  for  the 
investigation  of  regional  events  at  wavelengths  shorter  than  their  characteristic  dimension  and 
to  provide  a  basis  for  source  retrieval  from  sparsely  recorded  historic  events.  Thus,  we  will 
begin  with  a  brief  review  of  the  types  of  seismograms  available  at  the  Seismological  Laboratoryat 
Caltech  for  the  last  50  years  of  recording  and  how  they  could  be  used  to  full  advantage. 

Presently,  there  are  very  few  on-scale  recordings  of  major  earthquakes  (magnitude  >  5.8)  at 
near-regional  ranges  of  10  to  500  km  except  for  a  small  number  of  low-gain  (gain  •  100)  Wood- 
Anderson  (WA)  type  seismograms.  Moderate  earthquakes  in  southern  California  with  magnitudes 
between  4  to  5.5  are  usually  recorded  on  one  or  more  of  the  following  instruments  in  Pasadena: 

1 .  The  norma!  WA  instrument  (r0  »0.8  sec.)  with  a  gain  of  2800.  WA  seismograms  provide  the 
basis  for  the  Richter  scale  and  the  local  magnitude,  ML\  2.  The  long-period  (T#  »  6  sec.)  torsion 
instrument  which  operates  at  a  gain  near  900;  3.  The  Press-Ewing  (rc  -  30  sec.)  instrument 
(LP  30-90)  which  operates  at  a  gain  of  about  2250.  The  instrumental  responses  for  these  three 
systems  are  given  In  Figure  1 .  Two  examples  of  the  tangential  component  of  the  low-gain  WA  are 
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Figure  1.  The  Instrument  responses  of  the  Instruments  recording  the  events  used  In  this  study. 


50 


Figure  2.  Comparison  of  aftershocks  to  main  events  on  different  Wood-Anderson  (WA)  torsion 
mstnjments.  Top:  100X  WA  record  of  the  Coyote  Mountain  event  (filtered  through  a  0.5  sec.  up  and 
down  triangle)  and  an  aftershock  of  the  1954  Borrego  Mountain  earthquake  recorded  on  the  long-period 
WA  instrument.  Bottom:  100X  WA  record  of  the  San  Miguel  earthquake  and  one  of  Its  aftershocks 
recorded  on  a  high-gain  (gain  is  about  900)  long-period  WA  torsion  instrument.  (San  Miguel  records  were 
provided  by  Joanne  Yoshimura). 
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EVENT  LOCATIONS 
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Figure  3.  Map  showing  locations  ot  Pasadena.  PAS.  and  events  used  in  this  study. 
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AFTERSHOCKS  OF  IMPERIAL  VALLEY  1979  EQ 


Figure  4.  Three -components  seismograms  of  the  aftershock  sequence  recorded  on  the  long-period 
Press-Ewing  30-90  instrument  at  Pasadena.  Note  the  difference  In  waveforms  and  relative  anplitudes  of 
the  three  boxed  aftershocks  A,  B,  and  C  between  components.  The  event  locations  are  given  In  Table  t. 


53 


displayed  in  Figure  2  along  with  examples  of  aftershock  tangential  recordings  on  the  long-period 
(T0  =  6  sec.)  torsions. 

In  the  lower  comparison  we  display  a  magnitude  6.5  Northern  Baja  (Feb.  15,  1956)  earth¬ 
quake  as  recorded  on  the  100X  WA  and  one  of  its  aftershocks  (magnitude  >5.0)  on  the  high  gain 
long-period  torsion.  In  this  example  we  chose  an  aftershock  that  looks  similar  to  the  main  event. 
Note  that  the  main  event  is  probably  100  times  stronger  at  the  longer  periods  but  because  of  dif¬ 
ferences  in  instrumental  response  they  appear  similar.  Apparently,  a  complete  trade-off  between 
shifts  in  corner-frequency  of  the  source  spectrum  and  instrument  response  has  occurred.  The 
upper  two  records  are  not  so  similar  but  these  two  events  are  located  at  different  distances  from 
Pasadena.  The  bottom  event  is  actually  an  aftershock  from  an  earlier  event,  the  1954  Borrego 
Mountain  earthquake  and  the  top  event  is  the  1969  Coyote  Mountain  earthquake  recorded  on 
the  100X  WA  filtered  through  a  0.5  sec.  up  and  down  triangle.  Figure  3  shows  the  location 
of  all  events  used  in  this  study.  If  we  knew  the  mechanisms  of  these  aftershocks  we  probably 
could  make  some  good  estimates  of  the  faulting  parameters  of  the  main  events  by  using  the 
well-known  empirical  Green’s  function  approach  (Hartzell,  1978).  This  approach  has  proven  very 
useful  in  strong-motion  simulations  of  main  events  from  their  aftershocks  by  assuming  that  the 
aftershock  and  the  main  event  have  the  same  mechanisms.  The  strength  of  this  procedure  lies  in 
eliminating  the  path  effects  by  assuming  they  are  included  in  the  aftershock  record  automatically. 
The  biggest  short-coming  of  this  method  is  having  an  aftershock  with  the  same  orientation  as  the 
main  event  at  the  appropriate  depth  and  range. 

Generally  we  do  not  know  the  mechanisms  of  many  of  these  moderate  aftershocks  because 
they  are  too  large  for  local  arrays  and  too  small  for  global  networks.  In  addition,  waveforms  of 
different  aftershocks  are  often  so  different,  even  if  they  are  spatially  close  to  each  other,  that  it  is 
not  easy  to  determine  which  aftershocK  can  be  used  as  an  empirical  Green’s  function  to  simulate 
the  main  event.  Examples  of  this  catagory  of  events  are  the  series  of  aftershocks  of  the  1979 
Imperial  Valley  earthquakes  recorded  on  LP  30-90  at  Pasadena  (Figure  4).  We  would  like  to  call 
attention  to  the  three  aftershocks  labelled  A,  B  and  C  where  motions  for  all  three  events  are  about 
the  same  size  on  the  vertical  (UP)  component  but  distinctly  different  on  the  horizontals  (EW  and 
NS).  Differences  In  waveforms  are  also  apparent  on  all  three  components.  Since  the  paths  from 
the  Imperial  Valley  to  Pasadena  are  essentially  the  same  for  the  three  events  we  would  conclude 
that  the  source  characteristics  (depths  and  faulting  parameters  such  as  strike,  dip  and  rake)  must 
be  different.  As  we  will  discuss  shortly,  event  A  is  probably  a  mid-depth  normal  dip-slip  event, 
whereas  event  B  is  a  shallow  normal  dip-slip  event.  Event  C  is  a  deep  strike-slip  event. 

In  order  to  study  these  seismograms  in  detail,  we  digitized  and  rotated  the  NS  and  EW  com¬ 
ponents  to  obtain  the  tangential  (SH  waves)  and  vertical-radial  (P-SV  waves)  as  displayed  in 
Figure  5.  At  these  periods  the  motion  appeal’s  to  be  well-behaved  in  that  the  P  -waves  are  not 
apparent  on  the  tangential  component.  Particle  motion  studies  of  the  type  discussed  by  Vidale 
(1986)  conducted  on  these  recordings  Indicate  that  the  first  10  sec.  of  record  is  consistent  with 
P  -waves  and  diffracted  sv  -waves  followed  by  Rayleigh  motion.  Similar  analysis  of  filtered 
rotated  torsion  records  indicates  that  separation  of  the  P-SV  and  SH  system  oocurs  down  to 
periods  of  about  1  second.  Thus,  it  appears  that  two-dimensional  models  may  prove  effective  in 
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Figure  5.  Rotated  seismograms  of  the  three  aftershocks  shown  In  Figure  4.  Maximum  amplitudes 
are  given  in  mm. 
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VARIATIONS  OF  GREEN’S  FUNCTIONS 
WITH  DISTANCE 


Figure  6.  Preferred  velocity  model  end  the  corresponding  Green 's  functions  gene  reted  with  e  source 
depth  of  7  km  at  the  location  of  the  1976  Brawley  earthquake  (P  -  262  km).  The  stations  ere  located 
at  SO  km  interval  with  the  last  station  as  the  Pasadena  station.  Notice  the  rapid  development  of  surface 
waves  as  soon  as  seismic  waves  leave  the  basin.  Maximum  amplitudes  are  given  In  mm  x  10"  ^ 
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Fig ura  7.  Comparison  of  the  profile  of  Green's  functions  corresponding  to  a  pure  strike-slip  source 
as  recorded  on  the  three  Instruments  displayed  earlier  In  Figure  1.  Amplitudes  are  given  in  mm  for  a 
moment  of  <*p  x  io25  dyne-cm. 
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removing  the  propagations!  distortions  so  that  source  retrieval  Is  possible.  For  example,  Helm- 
berger  and  Vidafe  (1987)  have  modeled  successful  propagational  features  of  2-D  structures  at 
short  distances. 

In  this  study,  we  will  derive  a  2-D  model  appropriate  for  a  cross-section  from  Imperial  Valley 
to  Pasadena  using  master  events,  known  sources,  and  the  forward  modeling  approach.  The 
resulting  Green's  functions  are  then  used  to  examine  some  of  the  more  important  events  located 
in  the  section  (Figure  3).  We  only  considered  the  tangential  motions  in  this  pilot  study  because 
less  computational  effort  is  involved  compared  to  solving  a  complete  P-S  V  system  recorded  on 
the  radial  and  vertical  components. 

Crustal  Cross-section  and  Green’s  functions 

Searching  for  suitable  Green’s  functions  by  trial-and-error  testing  can  be  a  time-consuming 
endeavor  but  the  basic  approach  has  proven  effective  in  previous  studies  (Vidale  et  al.  ,1985; 
Helmberger  and  Vidale,  1987).  Two  types  of  codes  were  employed,  namely  the  generalized  ray 
method  (GRT)  for  laterally  varying  layers  (Helmberger  et  al,  1985)  and  a  modified  finite-difference 
(FD)  technique  (Vidale  et  al,  1985).  The  first  method  is  analytical  and  can  be  used  effectively 
to  adjust  deeper  smoothly  varying  structures  for  proper  timing  and  critical  angle  positions.  The 
truncation  of  basins,  however,  requires  the  more  powerful  numerical  approach.  Some  of  the 
details  of  forward  modeling  straLg^s  for  this  situation  will  be  discussed  later  but,  first,  we  present 
in  Figure  6  our  best  fitting  model  and  corresponding  Green's  functions.  The  source  is  situated 
at  mid-depth  (7  km)  in  the  Imperial  Valley  basin  with  responses  given  at  50  km  intervals.  Note 
the  soft  basin  structure.  This  structure  is  responsible  for  the  rapid  development  of  surface  waves 
which  occurs  within  the  first  50  km  of  propagation.  The  waveform  appears  to  change  relatively 
slowly  in  traversing  the  remaining  path.  Note  the  ringing  waveforms  for  the  dip-slip  case  versus 
the  simpler  larger  amplitude  pulse  for  the  strike-slip  orientation  at  the  nearest  range.  This  effect 
can  be  explained  easily  by  the  corresponding  vertical  radiation  patterns  (Vidale  et  al.,  1985). 

Figure  7  shows  the  corresponding  synthetics  for  the  strike-slip  case  convolved  with  the  in¬ 
strument  responses  displayed  in  Figure  1.  Amplitudes  are  given  in  mm  for  a  moment  of  4*p  x 
to23  dyne-cm  where  p  is  the  density  at  the  source  location. 

Changing  the  source  positions  relative  to  the  basin  has  the  most  dramatic  effect  as  displayed 
in  Figure  8.  The  relatively  hard-rock  path  to  Pasadena  produces  little  waveform  distortions. 
Shooting  back  into  the  basin  produces  a  wave-train  similar  to  the  motions  observed  at  El  Centro 
for  Borrego  Mountain  events  (Helmberger  and  Vidale,  1987).  The  quick  change  in  waveforms 
observed  at  Pasadena  for  small  changes  in  source  positions  near  the  basin  boundary  is  the 
primary  reason  for  abandoning  a  flat-layered  structure. 

Long-Period  Modeling 

As  in  all  forward  modeling  attempts,  one  starts  with  the  best  geophysical  data  available 
for  constraining  the  initial  model.  Fortunately,  considerable  studies  have  been  conducted  in 
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Figure  8.  Profile  of  Green's  functions  with  the  source  now  outside  the  basin.  Little  waveform  distor¬ 
tion  is  observed  along  the  path  towards  Pasadena  (upper  four  traces)  while  very  complicated  waveforms 
develop  rapidly  in  the  basin  < lower  two  traces),  Indicating  the  important  effect  of  the  basin  edge  on  wave 
propagation  along  this  profile. 
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Figure  10.  Sensitivity  of  Green's  functions  to  depth.  Sources  are  put  in  lour  different  depths  at  the 
same  epicentre!  location  as  the  1976  Brawtey  event:  3.5,  7.0,  10.5  and  14.0  km.  In  general,  when  the 
source  is  still  in  the  basin,  more  surface  waves  are  generated  and  with  the  source  below  the  soft  and  slow 
sediments  in  the  basin,  the  waveforms  are  simpler  and  energy  drops  off  rapidly  with  time. 
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this  region.  For  the  Imperial  Valley  velocity  profile,  we  used  the  model  proposed  by  Fuis  et 
44.(1982)  and  McMechan  and  Mooney  (1980).  Just  outside  the  basin  we  used  the  results  from 
Hamilton  (1970)  who  investigated  the  Borrego  Mountain  aftershocks  with  controlled  calibration 
shots.  Hamilton's  results  suggest  a  thick  crust-mantle  transition  zone.  At  Pasadena  we  adopted 
the  model  proposed  by  Hadiey  and  Kanamori  (1977,  1979)  with  a  thick  crust  and  relatively 
sharp  Moho  transition.  These  vertical  profiles  were  assembled  by  allowing  for  a  gentle  dipping 
connecting  structure  and  produce  an  initial  cross-section  from  Imperial  Valley  to  Pasadena. 

Synthetics  for  this  starting  model  were  generated  and  compared  with  observed  records  for  the 
various  calibration  events  displayed  in  Figure  3.  For  example,  in  Figure  9  we  give  the  comparison 
between  predictions  for  the  Brawley  event  and  its  aftershock.  The  source  parameters  for  this 
event  were  determined  earlier  by  Heaton  and  Helmberger  (1978),  and  treated  as  known.  This 
assumption  allows  us  to  perturb  the  various  model  parameters  to  bring  about  agreements  in 
waveform  and  absolute  timing.  This  procedure  goes  relatively  fast  for  long-period  modeling  but 
becomes  increasingly  tedious  at  higher  frequencies. 

The  best  fitting  model  to  date  has  a  slow  mantle  velocity  of  4.28  km/sec,  a  northwest-thinning 
Moho-crustal  transition  layer  of  4.18  km/sec,  a  dipping  crustal-Moho  transition  layer  of  4.05 
km/sec,  a  dipping  lower  crust  of  3.78  km/sec  that  thins  out  to  the  northwest,  and  an  upper 
crust  of  3.38  km/sec  that  also  dips  northwest.  The  idealized  Imperial  Valley  basin  surface  has 
two  layers  of  very  slow  shear  velocities  of  1.0  krrVsec  and  2.34  km/sec,  corresponding  to  what 
Fuis  et  4/.  (1982)  described  in  their  P  -wave  refraction  profiles.  A  thin  layer  of  3  km/sec  that  thins 
out  at  the  edge  of  the  basin  lies  underneath  the  slow  sediments.  This  model  is  able  to  predict 
the  arrival  time  and  the  waveform  of  the  first  40  seconds.  The  secondary  arrival  occuring  about 
60  sec.  after  the  initial  motion  is  present  on  many  records  and  will  be  discussed  later. 

The  moment  required  to  match  the  amplitudes  is  2.1  x  1023  dyne-cm  which  can  be  com¬ 
pared  to  3.2  x  1023  dyne-cm  found  by  Heaton  and  Helmberger  (1978)  using  local  strong-motion, 
intermediate-period  records.  This  seems  like  a  reasonable  difference  since  lack  of  correction  for 
Q  and  scattered  energy,  instrumental  gain,  small  changes  in  strike  etc.  could  aocount  for  such 
a  difference. 

Finally,  with  the  model  adjusted  to  predict  synthetics  matching  the  master  event,  we  can 
easily  perform  depth  sensitivity  studies  as  displayed  in  Figure  10.  The  synthetics  are  generated 
with  the  source  at  a  distance  of  262  km  from  Pasadena,  so  the  source  is  inside  the  basin  for 
depths  of  3.5  and  7  km  and  below  the  basin  for  depths  of  10.5  and  14  km.  In  both  strike-slip 
and  dip-slip  mechanisms,  absolute  amplitudes  decrease  with  depth.  In  addition,  pure  strike-slip 
mechanism  results  in  higher  amplitudes  compared  to  a  pure  dip-slip  mechanism,  almost  twice 
the  amplitudes  at  all  four  depths.  In  general,  we  also  expect  to  see  a  shallow  source  exciting 
more  surface  waves  in  both  pure  strike-slip  and  pure  dip-slip,  and  a  deeper  source  tends  to 
show  less  complexity.  These  results  are  similar  to  those  found  in  fiat  layered  models.  After 
reviewing  mechanisms  for  events  in  this  region  from  Fuis  et  al.  1982;  Heaton  and  Helmberger, 
197d;  Johnson  and  Hadley,  1976;  Johnson  and  Hutton,  1982;  Liu  and  Helmberger,  1985,  a  total 
of  21  possible  orientations  for  each  event  was  considered.  The  best  fitting  combinations  are 
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Figure  11.  Long-period  modeling  of  nine  unknown  aftershocks  and  small  events  In  the  Imperial 
Valley.  Data  shown  are  recorded  on  the  LP  m  80-90  instrument.  The  modeled  mechanisms  are  also 
shown.  Results  indicate  that  event  "C  is  a  deep  strike-slip  event,  event  'B‘  a  shallow  dip-slip  event  and 
event  ‘A’  a  mid-depth  dip- slip  event. 


Figure  12.  Sensitivity  of  Green 's  functions  to  depth  with  the  source  at  Anza  distance.  The  waveforms 
are  simple  compared  to  those  displayed  in  Figure  10.  There  is  not  as  much  a  difference  In  the  complications 
of  waveforms  with  depth  as  in  the  previous  case. 
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Figure  13.  Long-period  modeling  of  three  Anza  events  and  two  aftershocks  of  the  1968  Borrego 
Mountain  event  and  one  of  the  1969  Coyote  Mountain  aftershock.  The  corresponding  focal  mechanism  is 
also  shown. 


displayed  in  Figure  11  with  depths,  magnitudes  and  moment  estimates  given,  see  Table  1  and 
2.  Note  that  event  C  as  discussed  earlier  is  modeled  as  a  deep  strike-slip  while  event  B  appears 
to  be  a  shallow  dip-slip. 

For  small  events  we  expect  source  durations  to  be  short  compared  to  the  30*90  instrumen¬ 
tal  response,  see  Figure  1,  and,  therefore,  we  should  be  able  to  fit  any  observed  long-period 
seismograms  by  a  linear  combination  of  the  Green's  functions  displayed  in  Figure  10.  A  total  of 
nine  events  with  unknown  source  parameters,  including  the  three  aftershocks  described  earlier 
and  the  Brawley  aftershock,  were  collected  as  a  data  set  to  which  the  Green's  functions  were 
applied.  These  events  are  listed  in  Table  1  where  the  locations  and  depths  are  given  according 
to  the  Caltech  catalog. 

The  sensitivity  of  Green’s  functions  to  depth  for  sources  at  Anza  appears  totally  different  from 
that  in  Imperial  Valley.  In  general,  the  waveform  is  very  simple  for  both  strike-slip  and  dip-slip 
mechanisms,  and  it  changes  only  slightly  with  depth  (Figure  12).  It  appears  that  the  deeper  the 
event,  the  simpler  the  waveform  it  created.  This  feature  proves  useful  in  fixing  the  depths  of 
events  as  did  the  amount  of  surface  wave  excitation  for  events  in  the  Imperial  Valley. 

In  order  to  test  the  derived  crustal  model  for  its  accuracy  outside  the  Imperial  Valley  sed¬ 
imentary  basin,  we  applied  the  same  process  to  three  events  at  Anza.  an  aftershock  of  the 
Coyote  Mountain  1969  event  and  two  aftershocks  of  the  1968  Borrego  Mountain  event,  see  Fig¬ 
ure  13.  The  mechanisms  used  to  model  the  Anza  events  were  again  collected  from  various 
studies  (Sanders  and  Kanamori,  1984;  Given,  1983)  and  adjustments  made  to  fit  the  data.  The 
fits  between  synthetics  and  data  are  reasonably  good  with  results  given  in  Table  3. 

Most  of  the  aftershocks  in  the  Borrego  region  preferred  mechanisms  similar  to  their  main 
events  with  results  given  in  Table  4.  The  mechanisms  for  the  main  events  are  given  by  (Thatcher 
and  Hamilton,  1973;  Ebel  and  HeJmberger,  1982;  Burdick  and  Mellman,  1976  and  Sanders  et  «J. 
,  1986).  The  waveforms  are  relatively  simple,  as  shown  In  Figure  13,  compared  to  the  Imperial 
Valley  events  in  Figure  1 1 ,  suggesting  the  complications  of  waveforms  due  to  the  existence  of 
the  basin. 

Strong  Motions  from  Imperial  Valley  Events 

One  of  the  primary  objectives  of  this  study  was  to  provide  Green’s  functions  appropriate 
for  strong  motions  in  the  frequency  domain  of  10  to  0.1  Hz.  Unfortunately,  we  do  not  have 
broad  band  data  at  these  intensity  levels  so  that  we  must  rely  on  recordings  from  the  100X  WA. 

At  10  Hz  the  problem  becomes  extremely  complex  and  the  motions  no  longer  separate  into 
the  P-SV  and  SH  systems.  Secondly,  we  no  longer  expect  events  of  this  magnitude  to  be 
as  simple  as  events  discussed  earlier.  Thus,  it  is  difficult  to  assess  the  adequacy  of  our  results 
since  we  can  only  compare  predicted  motions  at  Pasadena  based  on  independently  determined 
source  studies  from  near-in  data  at  Imperial  Valley.  Three  component  data  may  help  resolve  the 
source  properties  by  providing  more  data,  and  this  subject  will  be  addressed  in  a  later  effort.  At 
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this  stage,  we  will  limit  ourselves  to  the  1979  Imperial  Valley  earthquake  to  serve  as  an  example. 
The  secondary  energy  arriving  about  40  secs  after  the  Initial  motion  which  accompanies  many  of 
the  shallow  events  in  the  Valley  will  be  discussed  later. 

1979  Imperial  Valley  Earthquake 

Several  inversions  were  done  on  the  1 979  Imperial  Valley  mainshock.  The  general  consensus 
of  the  rupture  includes  an  initial  10  km  deep  epicenter  that  rupture  northwestward  along  the 
Imperial  Fault  at  a  rate  of  75%  of  the  shear  velocity,  then  rupture  continued  on  at  a  shallower 
depth  (Hartzell  and  Helmberger,  1982;  Hartzell  and  Heaton,  1983;  Olson  and  Aspel,  1982). 
Archuleta  (1984)  holds  a  slightly  different  conclusion  on  the  rupture  process,  with  an  initial  strike- 
slip  source  at  about  8  km  depth  rupturing  northwestward.  Subsequent  rupture  occurred  at  two 
main  locations,  both  at  about  10  km  depth,  with  a  minor  dip-slip  rupture  at  about  30  km  from  the 
epicenter  along  the  strike  of  the  fault  The  Hartzell  and  Helmberger's  model  (model  HH1)  and 
the  Hartzell  and  Heaton’s  model  (model  HH2)  are  very  similar,  while  Olson  and  Aspel’s  model 
(model  OA)  is  a  more  continuous  model,  which  can  be  simulated  using  nine  segments  of  rupture. 
We  attempt  to  model  this  1979  mainshock  by  treating  each  rupture  segment  in  each  inversion 
model  as  an  earthquake  source.  Using  the  same  mechanism  as  inverted  by  the  above  workers 
in  each  segment,  we  combine  our  synthetics  with  the  mechanisms  and  then  add  the  segments 
up  according  to  the  corresponding  delay  time  along  the  fault.  Simulations  of  such  strong  motions 
appropriate  for  the  Pasadena  torsion  are  shown  in  Figure  14  with  the  corresponding  inversion 
models. 

A  comparison  of  the  filtered  data  with  the  synthetic  predictions  is  displayed  In  Figure  15. 
The  synthetic  responses  shown  in  Figure  14  require  a  time-derivative  to  compare  with  the  corre¬ 
sponding  WA  record  in  displacement  which  tends  to  emphasize  the  high  frequency  tails  such  as 
in  the  Olsen-Apsel  model.  In  general,  all  of  these  models  display  some  merit  although  it  appears 
that  the  two  models  on  the  left  fit  the  waveform  data  somewhat  better.  The  important  point  in  this 
comparison  is  not  which  model  fits  better  but  that  complex  earthquakes  (multiple-ruptures)  can 
be  distinguished  from  simple  events.  Thus,  this  type  of  regional  data  from  historic  events  can  be 
used  to  help  delineate  rupture  patterns  along  important  fault  segments. 

Discussion  and  Conclusion 

Many  of  the  events  in  the  Imperial  Valley  data  set  have  a  secondary  arrival  that  is  about  60 
sec  late  (Figure  11),  which  we  have  neglected.  There  are  several  possibilities  for  the  secondary 
arrivals:  (1).  Source  structural  effects  which  includes  three-dimensional  scattering  due  to  the 
Imperial  Valley  basin  structure,  (2).  Source  parameters  effects  which  suggests  double  events  or 
complicated  faulting  mechanisms,  (3).  Receiver  structural  effects  which  are  local  effects  due  to 
receiver  being  in  or  near  a  basin,  (4).  Path  complications  which  suggest  structural  effects  along 
the  same  path  that  are  not  in  the  present  model. 

Source  structural  effects  would  result  In  waves  arriving  at  the  receiver  along  different  az¬ 
imuths.  We  determined  by  complex  polarization  studies  (Vidale,  1986)  on  the  three  components. 
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Figure  15.  Corrparison  of  4  sec  filtered  tangential  motion  against  predicted  synthetics. 
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that  the  late  arrival  was  traveling  along  approximately  the  same  azimuth  as  the  main  arrival.  So 
we  believe  that  this  late  arrival  comes  mainly  from  the  same  ray  azimuth  as  the  main  arrival 
though  there  is  a  clockwise  rotation  after  the  first  arrival  at  Pasadena.  This  effectively  rules  out 
possibility  #1. 

The  three  events  we  studied  at  Anza,  the  Borrego  Mountain  distance  and  the  Coyote  Moun¬ 
tain  events  show  no  secondary  arrival  at  all.  The  arrivals  have  relatively  simple  waveforms.  The 
hypothesis  that  the  secondary  arrivals  on  the  Imperial  Valley  events  are  effects  of  the  local  re¬ 
ceiver  structure  is  then  ruled  out  because  of  the  absence  of  such  arrival  on  the  other  records  of 
events  outside  the  basin. 

We  are  now  left  with  possibilities  #2  and  #4.  It  is  fairly  unreasonable  to  attribute  double 
mechanisms  to  all  events  with  secondary  arrivals.  However,  shallow  events  have  secondary 
arrivals  while  deeper  events  do  not.  The  secondary  arrival  also  has  lower  frequency  than  the 
main  arrival  as  recorded  on  the  long-period  instrument  and  are  not  found  in  high  frequency 
records.  These  features  suggest  that  this  secondary  pulse  has  traveled  as  a  trapped  surface 
wave  in  the  surface  sediments  to  the  edge  of  the  basin  and  regenerates  into  a  normal  Love 
wave. 

A  number  of  numerical  profiles  with  sources  placed  at  various  depths  in  a  variety  of  basin 
models  were  generated  to  test  the  above  hypothesis.  When  the  source  is  situated  in  the  sed¬ 
iments  and  when  the  basin  ends  sharply  the  secondary  arrival  becomes  particularly  strong  as 
displayed  in  Figure  16.  Basin  models  with  gentle  dipping  edges  do  not  show  the  secondary  ar¬ 
rival  and  apparently  scatter  the  surface  waves  at  lower  ray  parameters,  probably  teleseismically 
as  found  in  the  study  by  Stead  and  Helmberger  (1988). 

The  broad-band  responses  displayed  In  Figure  16  do  not  contain  the  instrument  and  suggest 
that  the  secondary  arrival  is  not  depleted  in  high  frequency  as  observed.  Thus,  the  observed 
secondary  arrival  has  lost  its  high  frequency  by  attenuation  in  the  soft  sediments  or,  perhaps,  the 
source  excitation  has  very  low  stress  drop.  This  subject  is  best  pursued  with  broad-band  three 
component  array  data  and  will  be  addressed  in  a  subsequent  paper. 

In  conclusion,  we  have  demonstrated  the  effectiveness  of  using  events  with  known  source 
mechanisms  to  derive  a  regional  crustal  model  and  make  use  of  the  model  to  retrieve  source 
mechanisms  for  other  events.  This  method  is  proven  to  be  successful  in  predicting  SH  wave¬ 
forms  and  arrival  times  of  events  at  a  different  distance  from  the  receiver.  The  method  only 
requires  prior  knowledge  about  well  studied  events  and  the  availability  of  three-components  seis¬ 
mograms  with  well-recorded  events  of  unknown  sources.  These  requirements  can  be  met  easily 
since  we  have  abundant  stations  with  three-components  seismograms  that  have  recorded  many 
regional  events.  As  long  as  we  have  one  master  event  to  work  with  in  the  same  area,  source 
retrievals  and  strong  motion  simulations  or  predictions  can  be  done. 

We  have  also  demonstrated  the  sensitivity  of  various  Green’s  functions  to  depth  and  mech¬ 
anisms.  By  applying  the  depth  sensitivity  of  SH  waveforms,  we  can  now  have  an  estimate  of 
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the  depth  of  an  event  in  the  Imperial  Valley  recorded  at  regional  distance  by  the  rate  ol  drop¬ 
off  of  reverberations  in  the  record.  This  provides  a  fairly  powerful  tool  in  detail  modeling  of  an 
event  since  we  can  have  an  initial  estimate  of  the  depth  range  of  the  event.  With  such  a  de¬ 
tailed  model,  we  now  can  simulate  long  period  strong-motions  of  an  earthquake  with  any  source 
mechanism, location  and  magnitude  along  the  profile  from  Imperial  Valley  to  Pasadena. 
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Table  1.  List  of  earthquakes  used  in  the  Imperial  Valley  long  period  study 
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*  Brawley  Earthquake 


Table  2.  imperial  Valley  moments  obtained  by  fitting  long-period  and  high-frequency  SH 
and  the  corresponding  ml  and  depths  compared  to  the  catalogue  value. 

Mq  are  given  in  10JJ  dyne-cm.  Depths  are  given  In  km. 


Event 


K 

high  frequency 

long-period 

10.5 

15.3 

.6 


3.3 


4.8  4.8 


4.8 


5.2  5.0 


4.8  5.1 


.5  3.8 


Depth(km) 

Modeled 

catalog 

10.5 

6.14 

3.5 

4.27 

12.25 

5.08 
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Table  3.  Summary  of  source  mechanisms  fined  for 
three  Anza  events  in  the  same  magnitude  range  of  4  to  5. 


Event 

Depth 

Dip 

Rake 

Strike 

Date 

Time 

Longitude(W) 

760811 

15:24:55.42 

33®  28.0’ 

-116®  30.62’ 

14.0 

70® 

-86® 

45® 

780605 

16:03:03.72 

33®  25.21’ 

-116*  41.61’ 

14.0 

70* 

-00® 

152® 

700212 

04:48:42.26 

33®  27.21’ 

-116®  25.44’ 

81* 

D 

145® 

•Depths  are  given  in  km.  Strikes  are  measured  clockwise  from  North. 
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Table  4.  Summary  of  source  mechanisms  fitted  for 
aftershocks  of  the  1942,  54,  68  Borrego  Mountain 
and  1 969  Coyote  Mountain  earthquakes. 


Event 


Date 

Hrmn 

Sec 

Oct  20,1942 

1621 

57.0 

Oct  30,1942 

0535 

45.0 

Mar  23,1954 

0414 

50.0 

Apr  9,1968 

0800 

38.5 

Apr  9,1968 

1831 

38.5 

May  10,1069 

1440 

33.0 

[•a 


Dip 

Rake 

Strike 

1 

80° 

251® 

163" 

70® 

-86® 

56® 

80° 

251® 

163® 

45® 

00® 

0® 

80® 

251® 

163® 

80® 

251® 

163® 

•Depths  of  these  events  are  all  modeled  at  8  km,  except  for 
one  of  the  1942  aftershocks  (0535)  which  was 
modeled  at  12  km  depth.  Strikes  are  measured  clockwise 

from  North. 
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